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FOREWORD 


This  report,  "Prediction  of  Six-Degree-of-Freedom  Store  Separation 
Trajectories  at  Speeds  up  to  the  Critical  Speed,"  describes  a combined 
theoretical-experimental  program  directed  toward  developing  a computer 
program  for  predicting  the  trajectory  of  an  external  store  dropped  from 
an  aircraft  of  the  fighter-bomber  type  at  speeds  up  to  the  critical  speed. 
This  volume,  Volume  I - "Theoretical  Methods  and  Comparisons  with  Experiment," 
describes  the  theoretical  approach  and  presents  extensive  comparisons 
with  experimental  data.  The  second  volume,  Volume  II  - nUsers  Manual  for 
the  Computer  Programs,"  presents  detailed  instructions  on  the  use  of  the 
computer  programs. 

The  work  was  carried  out  by  Nielsen  Engineering  & Research,  Inc., 

850  Maude  Avenue,  Mountain  View,  California  94040,  under  Contract  No. 
F33615-71-C-1116.  The  contract  was  initiated  under  Project  8219,  Task 
821902,  of  the  Air  Force  Flight  Dynamics  Laboratory.  The  Air  Force  Project 
Engineer  on  the  contract  was  Mr.  Jerry  E.  Jenkins,  AFFDL/FGC.  The  report 
number  assigned  by  Nielsen  Engineering  & Research,  Inc.  is  NEAR  TR  37. 

The  authors  wish  to  thank  Mr.  Jenkins  for  his  assistance  during  the 
course  of  the  investigation.  Also,  they  would  like  to  thank  Mr.  Willard  E. 
Summers  and  Mr.  Robert  H.  Roberts  of  the  4T  Projects  Branch,  Propulsion 
Wind-Tunnel  Facility,  Arnold  Engineering  Development  Center,  for  the 
timely  performance  of  the  experimental  test  program. 
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ABSTRACT 


The  primary  objective  of  this  report  is  to  present  a prediction 
method  for  determining  the  forces,  moments,  and  six-degree-of-f reedom 
trajectories  associated  with  stores  released  from  single,  TER,  or  MER 
configurations  on  fighter-bomber  aircraft  up  to  the  critical  speed.  A sec- 
ondary objective  is  to  present  important  experimental  results  and  compari- 
son with  theory  from  an  extensive  wind-tunnel  test  program  designed  to 
provide  data  to  test  the  theory  through  systematic  measurements  of  flow 
fields,  store  load  distributions,  store  forces  and  moments,  and  captive- 
store  trajectories.  Extensive  comparisons  are  made  for  flow-field  angles, 
store  loading  distributions,  store  forces  and  moments,  and  store  six-degree- 
of-f reedom  trajectories.  Store  normal-force  distributions  and  side- 
force  distributions  are  presented  for  single  stores  as  influenced  by  air- 
plane angle  of  attack,  vertical  location  under  the  aircraft,  and  presence 
of  the  pylon.  Similar  results  are  presented  for  interference  effects  in 
TER  configurations.  Interference  between  the  two  clusters  in  a MER  con- 
figuration is  shown.  Empennage  effects  are  presented  and  compared  with 
theory  for  a range  of  roll  angles  and  angles  of  attack.  Comparisons  are 
presented  between  captive-store  trajectories,  and  those  predicted  by  the 
six-degree-of-f reedom  program.  Cases  are  shown  for  a single  store  with- 
out empennage  and  with  empennage.  A case  of  a store  released  from  a TER 
configuration  is  also  compared  with  the  captive  store  trajectory.  Generally, 
the  method  does  a good  job  of  predicting  detailed  store  loading  and 
six-degree-of-freedom  trajectories.  It  should  serve  as  a useful  preliminary 
design  tool. 
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velocity  components  of  store  center  of  moments 
in  inertial  system 

jt»  : 
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acceleration  components  of  store  center  of  moments 
in  inertial  system 

Poo 

free-stream  mass  density 

dihedral  angle  of  a wing  segment,  positive  as 
shown  in  figure  2 

A0 

potential  due  to  one  source  strip 

orientation  of  store  tail  fins  with  respect  to 
xs,ys,zs  coordinate  system  (fig.  10) 

<t> 

Pt 

pylon  thickness  envelope  slope  measured  positive 
as  shown  in  figure  6 

wing  thickness  envelope  slope  measure  positive  as 
shown  in  figure  5 

sweep  angle  of  a horseshoe  vortex  bound  leg 
measured  in  wing  or  pylon  chordal  plane,  positive 
swept  back  on  left  wing  panel  or  pylon 

sweep  angle  of  a source  strip  (figs.  5 and  6) 

ip 

p 

sweep  angle  of  a horseshoe  vortex  bound  leg 
measured  in  wing  or  pylon  planform  plane  positive 
swept  back  on  left  wing  panel  or  pylon 

y,e,«t 

store  yaw,  pitch  and  roll  angles  specifying  angular 
orientation  of  store  x,y,z  coordinate  system 
relative  to  £,rj,£  inertial  system,  see  figure  9 

XV 


A'P,A6,  A<I> 


changes  in  store  yaw,  pitch,  and  roll  angles  from 
carriage  position  values 


¥,e,<t> 


time  rates  of  change  of  ¥,0,  and  <t> 


xvi 


PREDICTION  OF  S IX -DEG RE E-OF- FREEDOM  STORE  SEPARATION 
TRAJECTORIES  AT  SPEEDS  UP  TO  THE  CRITICAL  SPEED 

Volume  I - Theoretical  Methods  and  Comparisons 
With  Experiment 


1.  INTRODUCTION 

This  report  is  the  final  technical  report  describing  a combined 
theoretical-experimental  program  which  has  been  carried  out  with  the 
objective  of  developing  a method  for  predicting  the  six-degree-of-f reedom 
trajectories  of  external  stores  separated  from  aircraft  of  the  fighter- 
bomber  type  at  speeds  up  to  the  critical  speed.  The  work  documented  in 
this  report  represents  an  extension  and  improvement  of  the  work  reported 
in  reference  1.  In  that  work  the  method  was  limited  to  three  degrees  of 
freedom.  The  pylon  was  not  accounted  for,  nor  was  the  TER  or  MER  rack. 

The  pylon  and  rack  are  accounted  for  in  the  present  work.  In  addition  a 
new  vortex-lattice  model  has  been  adopted  which  allows  for  mutual  inter- 
ference between  wing  and  pylon,  as  well  as  wing  dihedral,  two  things  which 
were  not  included  in  the  vortex-lattice  model  of  reference  1.  A new  wing- 
pylon  thickness  model  has  been  adopted.  The  six-degree-of-freedom  equations 
of  motion  used  in  the  trajectory  calculations  allow  for  mass  and  inertia 
asymmetries  of  the  store. 

The  experimental  program  carried  out  in  conjunction  with  the  present 
theoretical  program  was  directed  toward  supplementing,  the  data  obtained 
during  the  experimental  program  conducted  simultaneously  with  the  theo- 
retical work  of  reference  1.  These  earlier  data  are  tabulated  in 
reference  2 and  the  data  obtained  in  the  present  program  are  tabulated  in 
reference  3.  Some  of  the  data  are  summarized  in  graphical  form  in.  four 
AEDC  reports,  references  4 through  7.  The  overall  experimental  program 
was  designed  such  that  interference  effects  could  be  isolated  through  a 
component  by  component  build-up  of  the  parent  aircraft.  At  each  stage  of 
the  build-up  flow  field  survey,  pressure  distribution,  and  force  and 
moment  data  were  taken.  A representative  set  of  trajectories  was  also 
obtained  for  comparison  with  theory. 

The  next  sections  of  this  report  will  discuss  the  general  method  of 
approach  to  the  interference  problem,  the  mathematical  models  for  the 
aircraft  components,  the  calculation  of  the  flow  field,  the  calculation 
of  the  forces  and  moments,  and  finally  the  calculation  of  the  trajectory. 
Methods  which  have  not  changed  from  the  work  of  reference  1 will  be 
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summarized  and  new  methods  will  be  described  in  detail.  Among  the 
latter  are  the  new  wing-pylon  model,  the  empennage  force  and  moment  calcu- 
lation, and  the  equations  of  motion  with  mass  and  inertia  asymmetries. 

Finally,  extensive  comparisons  are  made  with  selected  results  from 
the  wind-tunnel  program  for  flow-field  angles,  store  loading  distributions, 
store  forces  and  moments,  and  store  six-degree-of-freedom  trajectories. 

Store  normal-force  distributions  and  side-force  distributions  are  presented 
for  single  stores  as  influenced  by  airplane  angle  of  attack,  vertical 
location  under  the  aircraft,  and  presence  of  the  pylon.  Similar  results 
are  presented  for  interference  effects  in  TER  configurations  for  both 
the  bottom  store  and  the  shoulder  stores.  In  addition,  the  effect  of 
the  TER  rack  is  shown.  Interference  between  the  two  clusters  in  a MER 
configuration  is  shown.  Empennage  effects  are  presented  and  compared  with 
theory  for  a range  of  roll  angles  and  angles  of  attack.  A special  assess- 
ment is  made  of  the  effects  of  store  boattail  on  store  loading  distributions, 
forces  and  moments,  and  boundary-layer  separation  characteristics.  Also, 
the  effects  of  store  boattail  on  empennage  forces  and  moments  are  investi- 
gated. Comparisons  are  presented  between  captive-store  trajectories,  and 
those  predicted  by  the  six-degree-of-freedom  program.  Cases  are  shown 
for  a single  store  without  an  empennage  and  with  an  empennage.  A case  of 
a store  released  from  a TER  configuration  is  also  compared  with  the 
captive-store  trajectory.  Using  the  computer  program,  the  effect  of 
damping  is  shown  for  a finned  store. 

2.  GENERAL  METHOD  OF  APPROACH 

The  determination  of  the  forces  and  moments  acting  on  an  external 
store  due  to  the  complete  aircraft  is  a complicated  problem  in  aero- 
dynamic interference.  The  approach  used  herein  is  to  make  a first  approxi- 
mation to  the  gross  interference  effect  termed  the  primary  interference » 
and  then  to  consider  the  small  residual  interference  as  an  additional 
interference.  The  primary  interference  can  be  viewed  as  the  first  term 
in  an  iterative  procedure.  Consider  the  entire  airframe  except  the 
store  in  question  which  produces  a nonuniform  flow  field  at  the  position 
to  be  occupied  by  the  store.  If  we  now  place  the  store  in  this  nonuniform 
flow  field  and  determine  its  reaction,  we  have  accounted  for  primary 
interference.  If  now  we  consider  the  effect  of  the  store  on  the  airframe, 
additional  velocities  are  induced  normal  to  wing,  fuselage,  and  pylon 
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which  must  be  cancelled  by  additional  source  and  vortex  distributions. 

This  interference  is  termed  additional  interference  and  depends  on  the 
store  location. 

In  the  present  method,  the  pylon  is  considered  part  of  the  airframe 
so  that  the  primary  interference  contains  most  of  the  interference  even 
for  the  attached  store.  The  success  of  this  method  is  in  proportion  to 
the  degree  to  which  the  primary  interference  encompasses  all  the  significant 
effects  so  that  the  additional  interference  can  be  neglected.  Methods 
for  calculating  additional  interference  have  been  developed  in  reference  1, 
where  it  is  shown  that  the  effect  of  the  additional  interference  on  the 
store  is  small.  However,  the  effect  of  the  additional  interference  on 
the  airframe  may  not  be  small. 

An  important  advantage  realized  by  making  the  additional  inter- 
ference effect  on  the  store  negligible  is  that  the  necessity  of  having 
to  solve  the  complete  mutual  interference  problem  for  the  store  and 
airframe  is  eliminated.  The  matrix  necessary  to  solve  the  complete  mutual 
interference  problem  requires  store  control  points  at  about  20  axial 
stations  and  at  a number  of  azimuthal  stations  in  addition  to  the  control 
points  on  the  airframe.  Elimination  of  those  additional  control  points 
probably  reduces  the  time  to  calculate  the  interference  field  to  about 
1/4  of  what  it  would  be  for  a complete  mutual  interference  calculation. 
Another  even  more  significant  fact  is  that  the  airframe  singularity  distri- 
butions yielding  the  primary  interference  field  do  not  change  with  store 
position  so  that  it  is  not  necessary  to  re-solve  the  interference  problem 
at  each  point  in  the  store  trajectory.  As  a consequence,  it  is  possible 
to  calculate  a store  trajectory  on  existing  computing  machines  in  a 
reasonable  time. 

In  the  method  described  in  this  report  the  following  aircraft  com- 
ponents are  accounted  for: 


(1) 

Fuselage  volume 

(2) 

Fuselage  angle  of  attack 

(3) 

Wing  thickness 

(4) 

Wing  angle  of  attack, 

camber,  and  twist 

(5) 

Pylon  thickness 

(6) 

Pylon  loading  due  to 

induced  sidewash  field 
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(7)  Rack  volume 

(8)  Store  volume 

Fuselage  volume  is  taken  into  account  by  approximating  the  shape  by  an 
axisymmetric  body  and  representing  this  body  by  a series  of  three- 
dimensional  point  sources  placed  on  the  body  axis.  The  volumes  of  the 
rack  and  all  of  the  stores  present  are  accounted  for  in  the  same  manner. 
Fuselage  angle  of  attack  is  taken  into  account  as  a simple  Beskin  upwash. 
Wing  and  pylon  thickness  are  accounted  for  by  representing  their  thickness 
distributions  by  a distribution  of  three-dimensional  source  panels.  The 
wing  angle  of  attack,  camber,  and  twist  and  the  pylon  loading  are  repre- 
sented by  a vortex  lattice. 

Certain  interference  effects  between  aircraft  components  are 
accounted  for  in  the  wing-pylon  vortex-lattice  calculation.  These  are 

(1)  Fuselage  in  uniform  flow  on  wing  and  pylon 

(2)  Stores  in  uniform  flow  on  wing  and  pylon 

(3)  Rack  in  uniform  flow  on  wing  and  pylon 

(4)  Wing  thickness  on  pylon 

(5)  Pylon  thickness  on  wing 

The  interference  velocities  induced  normal  to  the  wing  and  pylon  are 
treated  as  an  induced  camber  and  are  included  in  the  boundary  condition 
used  in  determining  the  strengths  of  the  vortices.  This  represents  a 
first  iteration  on  the  complete  mutual  interference  problem. 

A second  iteration  would  be  to  allow  the  wing  and  pylon  to  induce 
velocities  in  the  regions  of  the  fuselage,  rack,  and  stores  thus  causing 
them  to  be  in  a nonuniform  flow  field.  Since  the  source  distributions 
used  to  model  these  bodies  are  for  a uniform  flow,  higher  order  singu- 
larities would  have  to  be  used  to  account  for  this  nonuniform  flow.  The 
influence  of  these  higher-order  singularities  would  die  off  more  rapidly 
with  distance  and  thus  it  is  considered  a second-order  effect  and  is 
not  included. 

With  the  source  and  vorticity  distributions  which  represent  the 
various  aircraft  components  determined,  the  trajectory  can  be  calculated. 
This  requires  that  the  forces  and  moments  acting  on  the  ejected  store, 
including  damping,  be  calculated  at  each  point  in  the  trajectory.  This 
is  done  by  calculating  the  three-dimensional  velocity  field  in  which  the 
store  is  operating  as  seen  by  the  store.  To  do  this  the  store  is  removed 
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from  the  flow  field  and  the  perturbation  velocity  field  due  to  the  air- 
craft components  is  calculated  at  a number  of  points  in  the  region  the 
store  had  occupied.  The  store  is  then  placed  back  in  this  velocity  field. 
The  free-stream  velocity  the  store  sees  is  added  to  the  perturbation  field 
as  are  the  store  angular  velocities.  The  velocity  field  is  calculated  over 
the  store  length  in  this  manner.  The  axial  distribution  of  store  body 
forces  and  moments  is  calculated  from  this  field  and  then  integrated  over 
the  store  length.  A similar  procedure  is  used  to  determine  the  empennage 
forces  and  moments  with  the  velocity  field  calculated  at  a number  of  points 
over  each  tail  fin.  The  method  presently  treats  planar  and  cruciform 
empennages. 


3.  COMPRESSIBILITY  CORRECTION 

The  compressibility  correction  used  in  the  present  work  is  the  same 
as  that  used  in  the  work  of  reference  1 (see  section  3 of  that  report) . 
The  compressible  space  (x,  y,  z)  and  the  incompressible  space  (x'jy'jZ*) 
are  related  by  the  following  transformation 


x' 


z ' = z 


and  the  velocities  are  related  by 


u = 


(1) 


(2) 


Also  in  reference  1 it  is  shown  that  the  incompressible  angle  of  attack, 
a 1 , is  given  by 


a 1 = 3a  (3) 

The  procedure  used  to  calculate  the  perturbation  velocity  field  used 
in  determining  the  forces  and  moments  acting  on  a store  in  the  presence 
of  a parent  aircraft  will  now  be  briefly  described.  The  first  step  is  to 
determine  the  source  and  vorticity  distributions  which  model  the  various 
components  of  the  parent  aircraft.  To  do  this  the  entire  configuration 
(fuselage,  wing,  pylon,  rack  and  store(s))  is  transformed  to  an  equivalent 
one  for  incompressible  flow.  A common  coordinate  system  whose  origin  is 
at  the  fuselage  nose  is  adopted  and  all  components  are  located  in  this 
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system.  For  a given  Mach  number  the  transformation  given  by  equation  (1) 
is  applied  to  determine  the  equivalent  incompressible  configuration.  This 
transformation  results  in  an  increase  in  the  lengths  of  the  fuselage, 
rack,  and  store ( s) , a stretching  of  the  wing  and  pylon  chords  and  an 
increase  in  their  sweep  angles.  In  addition,  if  the  wing,  rack,  or 
store (s)  is  at  incidence  with  respect  to  the  fuselage  longitudinal  axis, 
the  incidence  angle  is  reduced  by  the  factor  3.  The  fuselage  angle  of 
attack  is  also  reduced  by  this  factor.  Source  distributions  are  then 
obtained  to  represent  the  volume  distributions  of  the  equivalent  fuselage, 
rack,  and  stores.  The  vortex-lattice  and  source-strip  distributions  are 
laid  out  on  the  equivalent  wing-pylon  combination  which  is  stretched  in  the 
chordwise  direction.  Using  the  common  fuselage  coordinate  system,  the  pertur- 
bation velocities  induced  by  the  equivalent  configuration  at  points  of 
interest  in  the  incompressible  space  are  calculated  by  considering  each 
component  separately  and  then  adding  up  all  their  effects.  These  incom- 
pressible perturbation  velocities  are  then  transformed  back  to  compressible 
space  using  equation  (2) . 

In  calculating  the  store  forces  and  moments,  these  velocities  are 
calculated  in  the  region  the  store  longitudinal  axis  and  tail  fins  would 
occupy  if  the  store  were  in  the  flow  field.  This  perturbation  velocity 
field  is  then  resolved  into  the  store  coordinate  system. 

4.  MATHEMATICAL  MODELS  FOR  DETERMINING  FLOW  FIELD 

The  methods  to  be  presented  for  determining  the  forces  and  moments 
acting  on  a store  require  a detailed  knowledge  of  the  flow  field  in  the 
region  of  the  store  in  its  absence.  In  order  to  calculate  this  flow 
field  all  of  the  components  of  the  parent  aircraft  must  be  modeled  mathe- 
matically. It  is  the  purpose  of  this  section  of  the  report  to  describe 
these  models.  First,  the  flow  model  for  the  fuselage  and  stores  is 
described.  Next  the  wing-pylon  model  is  presented  and  finally,  an 
approximate  model  for  a rack  is  discussed.  All  of  these  models  are  for 
the  equivalent  incompressible  configuration. 

4.1  Fuselage  and  Store  Flow  Model 

The  flow  model  used  to  represent  the  volume  distributions  of  the 
equivalent  incompressible  fuselage  and  stores  is  the  same  as  that  pre- 
sented in  section  4.1  of  reference  1.  The  shapes  are  approximated  by 
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axisymmetric  bodies  and  the  resulting  volume  distributions  represented 
by  a distribution  of  three-dimensional  point  sources  along  the  body 
longitudinal  axis.  The  source  strengths  are  determined  in  the  manner 
described  in  reference  1. 


With  the  source  strengths  determined,  the  perturbation  velocity 
field  due  to  the  body  can  be  calculated.  Consider  figure  1 which  shows 
a right-hand  x,  y,  z coordinate  system  fixed  in  the  body  nose  with 
positive  directions  as  shown.  Also  shown  are  the  axial  and  radial 
perturbation  velocities  at  a field  point,  u and  v^  respectively,  with 
positive  directions  as  shown.  From  reference  1 the  expressions  for 
these  velocities  produced  by  N sources  are 


where 


N 


u*  ( 


„ Q£(x*  - x*) 

x*’r*>  = L 1 — W7 

k-i  [<**  - + r*2l 


N 


Q*r* 


v*(x*,r*)  = Y -p r—^ 

k-i  [<**  - *t>s  * r*2]  7 


' ★ — 


= 


r , Jfc 


2 + z2 


u 


* = 


R 


u 

V 


V* 

r 


v 


V 


Q£  = 


47r^RV00 


(4) 


(5) 


and  where  £ is  a reference  length  which  for  convenience  will  be 
taken  as  the  length  of  the  incompressible  body.  The  quantities  Q* 

x . -i 

and  x£  are  the  source  strength  and  location,  respectively,  of  the  k 
point  source,  and  the  point  (x*,r*)  designates  the  field  point  in 
cylindrical  coordinates.  In  resolving  v^  into  components  in  the 
y and  z direction,  v and  w respectively,  proper  account  must  be  taken 
of  the  quadrant. 
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4.2  Wing-Pylon  Flow  Model 

In  order  to  represent  the  wing-pylon  combination  by  a potential  flow 
model,  two  types  of  singularities  are  used  simultaneously.  A three- 
dimensional  vortex  lattice  is  laid  out  on  the  surfaces  of  both  wing  and 
pylon  to  account  for  aerodynamic  loadings,  and  a distribution  of  a special 
type  of  sources  accounts  for  thickness  effects.  The  source  strengths  are 
proportional  to  the  thickness  envelope  slopes  of  the  wing  and  the  pylon. 

The  flow  tangency  boundary  condition,  including  all  wing-pylon  interference, 
is  applied  at  a finite  number  of  points  on  the  wing-pylon  combination 
resulting  in  a set  of  simultaneous  equations  from  which  the  vortex 
strengths  are  determined.  A Prandtl-Glauert  similarity  rule  (see  section  3) 
is  used  to  account  for  compressibility. 

4.2.1  Configuration  characteristics 

The  configuration  of  interest  is  a wing  with  one  pylon  under  each 
wing  panel.  The  right  wing  panel  with  pylon  is  the  image  of  the  left 
wing  panel  and  pylon.  The  various  wing-pylon  parameters  included  in 
the  computer  program  are  listed  below. 

Wing  Panels: 

Thickness  distribution:  Specified  at  large  number  of  chordwise 

locations  and  at  the  same  number  of  spanwise  locations  as  used  in  the 
vortex  lattice  which  represents  the  wing  as  a lifting  surface. 

Mean-camber  surface:  May  have  both  twist  and  camber. 

Leading-edge  shape:  Up  to  thirty  straight  line  segments  of  differing 

sweep. 

Trailing-edge  shape:  Same  as  for  leading  edge. 

Taper:  Controlled  by  shape  of  leading-  and  trailing-edges. 

Tips:  Straight. 

Dihedral:  Variable  across  wing  panel. 

Pylon: 

Thickness  distribution:  Same  method  of  description  as  for  wing 

panel . 

Mean-camber  surface:  Planar. 

Leading-edge  shape:  Straight  line  which  may  be  swept. 

Trailing-edge  shape:  Same  as  for  leading-edge  except  sweep  angle 

may  be  different. 
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Taper:  Uniform. 

Tip:  Parallel  to  pylon  root  chord  (pylon-wing  junction) . 

Orientation:  Vertical  and  streamwise. 

4.2.2  Vortex-lattice  model 

The  three-dimensional  horseshoe  vortex-lattice  approach  that  is  used 
to  represent  the  wing-pylon  lifting  surfaces  follows  the  methods  developed 
in  references  8 and  9.  In  what  follows,  the  theory  and  layout  of  the 
vortex  lattice  will  be  described  first.  The  flow  tangency  boundary 
conditions  specified  for  the  wing  and  pylon  will  show  how  mutual  inter- 
ference due  to  aerodynamic  loading  as  well  as  thickness  effects  are 
accounted  for.  The  theory  and  layout  of  the  sources  superimposed  on  the 
vortex  lattice  to  account  for  thickness  will  be  discussed  in  section  4.2.3. 

The  wing  vortex-lattice  method  used  in  the  previous  store  trajectory 
work,  reference  1,  has  been  extensively  modified  to  treat  a wing-pylon 
combination  as  an  integral  unit.  The  new  lattice  program  also  allows 
for  breaks  in  wing  sweep  and  dihedral.  Wing  twist  is  assumed  small  and 
is  accounted  for  only  in  the  flow  tangency  boundary  condition.  As  such, 
the  twist  is  not  included  in  the  geometry  of  the  wing-pylon  combination 
so  that  the  wing  tip  is  not  twisted  relative  to  the  wing  root  chord. 

The  winq  and  horseshoe  vortex  coordinate  systems  for  a wing  with  breaks 
in  leading  and  trailing  edge  sweep  and  dihedral  are  illustrated  in 
figure  2.  Figure  3 shows  the  vortex-lattice  arrangement  for  a swept 
wing  with  constant  dihedral  and  pylons.  The  wing  panels  and  pylons 
are  divided  into  trapezoidal  area  elements.  A horseshoe  vortex  is 
placed  in  each  area  element  such  that  the  spanwise  bound  leg  lies  along 
the  element  quarter  chord  and  its  trailing  legs  along  the  sides  of  the 
element.  The  trailing  legs  are  assumed  to  lie  in  the  plane  of  the  area 
element.  The  area  elements  in  each  chordwise  row  have  equal  chords  and 
spans.  In  the  spanwise  direction,  the  area  element  widths  need  not  be  equal 
to  allow  for  closer  spacing  where  large  spanwise  loading  gradients  exist. 

For  a wing  with  breaks  in  sweep,  the  area  elements  are  arranged  spanwise 
so  that  the  breaks  in  leading-  and/or  trailing-edge  sweep  lie  on  the  line 
formed  by  one  of  the  sides  of  a chordwise  row  of  elements.  A similar 
procedure  is  followed  when  there  are  breaks  in  dihedral.  The  wing-pylon 
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junction  is  also  made  to  lie  along  a common  boundary  between  two  adjacent 
rows  of  elements-  Consequently,  the  trailing  legs  of  the  upper  vortices 
on  the  pylon  lie  at  the  junction,  as  shown  in  figure  3.  The  length  of 
the  wing-pylon  junction,  the  pylon  root  chord,  need  not  be  equal  to  the 
wing  local  chord. 

On  the  wing,  the  flow  tangency  boundary  condition  should  be  expressed 
at  the  camber  surface  but  instead  it  is  applied  at  a finite  set  of  control 
points  given  by  the  midpoint  of  the  3/4-chord  line  of  each  area  element 
located  in  the  wing  chordal  plane.  This  is  the  planar  approximation. 

The  chordal  plane  of  wing  segment  i is  the  plane  containing  the  leading 
and  trailing  edges  of  the  wing  segment  and  making  dihedral  angle  0^  with 
the  Z = 0 plane.  On  the  camberless  pylon,  the  flow  tangency  condition 
is  applied  likewise  at  control  points  situated  in  its  chordal  plane 
determined  by  halving  the  thickness  envelope. 

The  boundary  condition  states  that  there  is  no  flow  through  the 
wing  and  pylon  surfaces  at  each  control  point.  The  condition  is  illus- 
trated in  figure  4 for  the  wing.  The  velocities  normal  to  the  wing  consist 
of  a component  of  the  free  stream,  perturbation  velocities  u,  v,  and  w 
induced  by  the  wing-pylon  horseshoe  vortex  system,  and  perturbation 
velocities  u^,  v^,  and  w^  induced  by  the  distribution  of  sources, 
accounting  for  thickness  superimposed  on  the  vortex-lattice.  Velocities 
ui,  v.,  and  w.  may  also  include  other  externally  induced  perturbation 
velocities  such  as  those  due  to  the  fuselage  and  stores.  With  M control 
points  on  the  left  wing  panel  and  MP  on  the  left  pylon,  the  boundary 
condition  on  the  left  wing  panel  is  given  by 


cos  0, 


v 


v . 


= (a  + a ^ ) cos  0^  + 


v 


(6) 
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the 


The  wing  angle  of  attack  is  a and  is  the  local  angle  of  attack 
due  to  wing  camber  and  twist.  Both  angles  are  assumed  small.  With 
camberless  pylon  at  zero  incidence,  the  pylon  boundary  condition  is 
written  for  MP  control  points  as: 


M 


M+MP 


- Y. 


n_ 

47 rV 


n=i 


v,n 


-i 


47 TV 


V 


n=M+i 


v,n 


v . 


V 


V=M+1,  M+2 , . . .M+MP 


(7) 


The  right-hand  side  of  equation  (6)  represents  the  free-stream  component 
and  the  externally  induced  perturbation  velocities  normal  to  the  wing 
chordal  plane.  The  right-hand  side  of  equation  (7)  consists  only  of  an 
external  perturbation  velocity  normal  to  the  pylon  chordal  plane  since 
there  is  no  pylon  incidence.  The  first  summation  on  the  left-hand  side 
of  the  equations  represents  the  perturbation  velocities  induced  by  vor- 
ticity  on  the  left  and  right  wing  panels.  The  second  summation  repre- 
sents the  velocities  induced  by  vorticity  on  the  left  and  right  pylons. 
The  functions  inside  the  summations  will  be  discussed  next. 

The  functions  F , F , and  F are  influence  functions  relating  the 
u v w 3 

perturbation  velocity  components,  induced  at  some  point  by  a horseshoe 
vortex,  to  its  circulation  and  the  coordinates  of  the  point  relative  to 
the  vortex.  The  relationship  is  obtained  from  the  Biot-Savart  law,  see 
reference  10.  For  example,  the  wash  velocities  at  a point,  induced  by 
a horseshoe  vortex  on  the  wing,  are 

u(x,y,z)  = Fu(x,y,z,s ,^,0) 


v (x,y, z) 


_T_ 

4 IT 


Fv(x,y,z,s,^,0) 


(8) 


w(x,y,z)  = Fw(x,y,z,s,0,0) 

The  influence  functions  depend  only  upon  the  coordinates  (x,y,z)  relative 
to  the  horseshoe  vortex  of  the  point  at  which  the  velocity  is  to  be  com- 
puted, and  the  vortex  span,  s,  bound  leg  sweep  angle  tf/y  and  wing  segment 
dihedral  angle,  0. 

Consider  first  a vortex  on  the  left  wing  panel.  The  coordinate 
system  is  shown  in  figure  2.  Since  the  horseshoe  vortex  lies  in  the 
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chordal  plane  (the  plane  of  the  wing  panel) , the  bound-leg  sweep  angle  ip 
used  is  the  angle  in  the  chordal  plane,  whereas  the  sweep  angle  is 
generally  specified  in  the  planform  plane,  ^ . The  two  are  related  by 

Jr 

the  dihedral  angle  as 


tan  ip  = 


tan  7p 

P 


COS  0 


(9) 


The  induced  velocity  is  computed  by  integrating  over  the  lengths 
of  the  two  trailing  legs  and  the  bound  leg.  For  a horseshoe  vortex  on 
the  left  wing  panel,  the  resulting  influence  functions  are  as  follows. 
The  backwash  influence  coefficient,  positive  forward,  is 


,(x,y,z,s  ,^,0) 

(z  cos  0 — v sin  0 

cos  f 

£x  cos  ip  - (y  cos  0 + z sin  0)  sin  7p 

2 + (z  cos  0 - y sin  0)2 

e 


i x + s tan  tI/)  sin  & + (y+s  cos  0)cos  it/  cos  0 + (z  + s sin  0)  cos  il/  sin  0 


£(x  + s tan  7p)  2 + (y  + s cos  0)  2 + (z  + s sin  0)  2 j 
(x  — s tan  tp)  sin  0 + (y  - s cos  0)  cos  tJ/  cos  0 + (z  - s sin  0)  cos  i/  sin  0 


|(x  - s tan  0)  2 + (y-s  cos  0)  2 + (z-s  sin  0) 2 j1/2 


} 


(10) 


The  sidewash  influence  coefficient,  positive  to  the  right,  is 
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Pv(x,y ,z,s,^,0) 


- z sin  il'  4 x cos  V/  sin  0 


[x  cos  if/  - (y  cos  0 4 z sin  0)  sin  0J2  + (z  cos  0 - y sin  0)2 

(x  + s tan  if/)  sin  # + (y  + s cos  0)  cos  jr  cos  0 4 (z  4 s sin  0)  cos  sin  0 


e 


1/2 


£(x  4 s tan  if/) 2 4 (y  4 s cos  0) 2 4 (z  4*  s sin  0)2 
(x  - s tan  il/)  sin  tf/ + (y  - s cos  0)cos  ^ cos  0+  (z  - s sin  0)  cos  if/  sin  0 


[(X  - s tan  ^)2  4 (y  - s cos  0)2  + (z  - s sin  0)2J 


1/2 


■} 


(z  - s sin  0) 


(y  - s cos  0 ) 2 4 (z  - s sin  0); 


. 1 - 


(x  - s tan  if/) 


£ (x  — s tan  if/)2  + (y-s  cos  0)  2 4 (z-s  sin  0)  2 J 


1 /aj 


(z  4-  s sin  0) 

(Y  4-  s cos  0)p  4 (z  4-  s sin  0) ; 


'1  - 


(x  4-  s tan  if/) 


[ (x  + s tan  if/)  2 4-  (y  4*  s cos  0)  2 4 (z  + s sin  0)  2J 


i/v 


(ii) 


The  downwash  influence  coefficient,  positive  downward,  is 


Fw(x,y,z,s,^,$) 


- x cos  if/  cos  0 4 y sin  if/ 


[x  cos  if/  - (y  cos  0 + z sin  0)  sin  if/ ]2  + (z  cos  0 - y sin  0) 


e 


i x 4 s tan  ^//)  sin  ^4  (y  + s cos  0)  cos  1!/  cos  0 4-  (z  4 s sin  0)  cos  if/  sin  0 


|(x  4 s tan  if/)2  4-  (y  4 cos  0)2  4-  (z  4 s sin  0)2] 


1/2 


(x  - s tan  if/)  sin  if/  + (v  - s cos  0)  cos  if/  cos  0 4 (z-s  sin  0)  cos  0 sin  0 
[(x  - s tan  if/)2  4 (y  - s cos  0)2  4 (z  - s sin  0)2J1//2 


} 


n - s cos  0) 


(y-s  cos  0 ) 2 4 (z  - s sin  0); 


(Continued  on  next  page) 
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r 

• <i 


( x — s tan  il/) 


£(x-s  tan  0)2  + (y  - s cos  0)2  + (z  - s sin  0)2]1/^?- 


(y  4-  s cos  0) 

(y  + s cos  0)2  + (z  + s sin  0) 2 


(x  4-  s tan  ^/) 

£(x  + s tan  0)2  + (y  + s cos  0)2  + (z  + s sin 


(12) 


For  the  case  of  unyawed  vortices,  0 = 0°,  the  influence  coefficients 
given  above  degenerate  to  the  equations  given  by  Blackwell,  reference  11. 

Once  the  influence  coefficient  for  the  vortex  on  the  left-wing  panel 
is  obtained,  the  coefficient  for  the  vortex  at  the  image  position  on  the 
right-wing  panel  can  be  obtained  by  changing  appropriate  signs  in 
equations  (10)  through  (12)  and  by  taking  into  account  the  change  in  the 
y coordinate-  Thus,  for  both  wing  panels, 

Fu  = Fu<x,yrz,s, *,«*>)  + Fu(x,yr,z,s,-^,-<J>)  (13) 


where  the  first  term  is  the  left  panel  vortex  contribution  and  the 

second  term  the  right  panel  contribution.  The  same  transformations 

hold  for  F and  F . 

v w 

In  similar  fashion,  the  contributions  of  the  vortices  on  the  pylons 
can  be  obtained  from  the  left  wing  panel  expressions.  Thus,  for  the 
pylons 


Fu  = Fu(x»y^»z»s.^>^  “ -tt/2)  + Fu  ( x , y r , z , s , —if/ , <(> , = v/2) 


(14) 


where  the  first  term  is  the  contribution  of  the  left  pylon,  which  has 
a "dihedral  angle"  of  -tt/2,  and  the  second  term  represents  the  right 
pylon  contribution.  The  same  transformations  again  hold  for  F^  and  F^. 

Equations  (6)  and  (7)  represent  a set  of  M+MP  simultaneous 

equations  in  which  the  unknowns  are  the  M+MP  values  of  circulation 

strength  T.  Therefore,  the  values  can  be  obtained  through  a matrix 

solution  for  a given  angle  of  attack  and  a specified  set  of  externally 

ui  vi  w 4 

induced  perturbation  velocities  — , — , ^ caused  by  wing  and  pylon  thick- 

oo  oo  00 

ness  effects  or  from  additional  external  sources  such  as  fuselage,  rack  and 
stores . 
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4.2.3  Thickness  model 

The  three-dimensional  source  method  used  to  account  for  wing  and 
pylon  thickness  is  known  as  the  spanwise  source  strip  method  and  is 
quite  similar  to  that  used  by  Fernandes  in  reference  12.  The  wing  will 
be  dealt  with  first.  The  coordinate  system  and  angle  convention  used  in 
connection  with  the  theory  and  layout  of  these  source  strips  is  shown 
in  figure  5 for  the  left  wing  panel  with  specified  thickness  envelope. 

It  is  noted  that  this  coordinate  system  is  the  same  as  the  one  used  for 
the  vortex-lattice  method,  figure  2.  The  incompressible  velocity  potential 
due  to  a surface  with  chord  | - Xg|  and  span  | - Y.J  shown  in 

figure  5 on  the  wing  may  be  obtained  from  reference  13  or  14  with  a 
simple  coordinate  transformation.  The  result  is 


A<£  _ _1_ 

V 2tt 
00 


tan  0t  dXx  dYx 

V (Xx  - X)2  + (Y  - Y]l  ) 2 + (z7  - Z)2 


(15) 


for  a surface  located  at  Z1  = constant  and  with  local  incidence  angle 
0^  measured  from  the  positive  X direction,  see  figure  5.  The  inte- 
grations are  performed  over  the  span  and  chord  of  the  surface  and 
coordinates  X,Y,Z  locate  the  field  point  at  which  the  velocity  potential 
is  determined.  To  simplify  the  integrations,  spanwise  strips  are  arranged 
on  the  wing  such  that  inclination  angle  0t  is  approximately  constant. 
These  strips  are  swept  with  the  wing  and  have  much  smaller  chord  than 
span  dimension.  Then  the  Xx  coordinates  of  points  in  the  narrow  strip 
are  approximated  by 

Xx  = Xc  + Yx  tan  Tu  (16) 

where  is  the  sweep  angle  measured  from  a direction  parallel  to  the 

positive  Y-axis,  see  figure  5,  and  Xc  is  the  value  of  X2  when  Yx 
is  zero. 

The  integral  then  simplifies  to 


M = 

V 


tan  <£tAXi 


2? r 


X 

/ 


dY, 


V (Xc  + Yx  tan  V - X)2  + (Y  - y[ ) 2 + (Zx  - Z) 


(17) 
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where  AX^  = Xg  - X^.  Note  that  Xg  and  X^  as  shown  in  figure  5 are 
negative.  The  integration  can  be  performed  and  noting  that  X.  = 180°  + ip. 
and  tan  <f>t  = -tan  0 ~ -0,  the  velocity  potential  due  to  one  strip  with 
span  | Y.  - Yg  | at  Z1  equals 


M _ 


0AX 


Vco  ~ 2tt 


— COS  Tp^  {^n  [v< 


(Xc  + Ya  tan  - X)2  + (Y  - Yg)2  + (Zx  - Z) 


cos  1p^ 


- + (X  - Xc)  sin  ip^  + Y cos  ip^ 


- ■*"  [ V(Xc  * Yb  tan  - XIs  4-  (Y  - Yb)  * + (Z,  - Z)  * - ^ - - 


(X  - X ) sin  ^ + Y cos  ^ 


} 


(18) 


where  on  the  left  wing  panel  Yg  and  Y^  are  negative.  Differentiating 
with  respect  to  X,Y,Z  to  obtain  the  perturbation  velocities  due  to  one 
strip  results  in  the  following  expressions: 


4u 

v«> 

II 

M.) 
v J 

0AXi 

= ^ cos 

.f 

"(xc  - 

■ X 

+ Y tan  ip, 

a X 

l 

AB 

Av 

V 

oo 

-4% 
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M\ 

VJ 

0AX± 

= 2tt  COS 

AW 

a / 

0AXi 

“ az  ( 

t) 

■ 2ir  cos 

CD 


!^)  + D sin 


B cos  ip.  Y - Y, 

r X D 


AB 


- Z)  -(Z,  - 


+ D COS 
CD  j 


AB 


i - zPl 

CD  J 


(19) 


where 


16 


A 


= -\fixl  - X + Y~  tan  0 . ) 2 + (Y  - YJ2  + (Z,  - Z)2 

1 C al  a 


COS  1p^ 


- + (X  - X ) sin  ip . -f  Y cos  tf/. 


= V(X^  - X + Ya  tan  0.)2  + (Y  - Y ) 2 + (Zx  - Z) 

f C ax  a 


- V< 


Xc  - X + Yb  tan  fL) 2 + (Y  - Yb) 2 + (Zx  - Z); 


(20) 


COS  Tp^ 


+ (X  - Xc)  sin  0^  + Y cos 


D = y(Xc  - X + Yb  tan  ^i>2  + (Y  - Yb>2  + (Zj  - Z)2 


Similar  expressions  can  be  derived  for  the  corresponding  source  strip  on 
the  right  wing  panel.  The  resulting  expressions  are  identical  to 
equations  (19)  and  (20)  except  that  Yq  and  are  interchanged,  since 

Ya  is  still  the  inboard  side  of  the  strip,  and  tan  ip  ^ is  replaced 
by  -tan  ip ^ and  cos  ip ^ is  replaced  by  -cos  ip On  the  right  wing 
panel  Yq  and  Y^  are  positive.  At  a field  point  (X,Y,Z)  the  pertur- 
bation velocity  components  —9  — induced  by  wing  thickness  are 

then  given  by  equations  (19)  and  ^20)  °°summed  over  all  wing  source  strips 
on  both  the  left  wing  panel  and  right  wing  panel. 

The  coordinate  system  and  angle  convention  used  with  the  pylon 
thickness  model  is  shown  in  figure  6.  This  coordinate  system  is  also 
the  same  as  the  vortex-lattice  coordinate  system.  Analogous  to  the 
wing  thickness  model,  the  incompressible  potential  due  to  a surface  at 
Yp  = constant  and  with  chord  | X^  - XQ | and  span  | - ZQ)  may  also 

be  obtained  from  reference  13  or  14  with  a simple  coordinate  transfor- 
mation. With  the  local  incidence  angle  measured  from  the  positive 

X-direction,  the  result  is 


M = 

~ 2tt 

00 


tan  0 dZ,  dX, 
Pt 


- X) 2 '+  (Yx  - Y)2  + (z[  - Z)2 


(21) 
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The  integrations  are  performed  over  the  span  and  chord  of  the  surface 
and  coordinates  X,Y,Z  locate  the  field  point  at  which  the  velocity 
potential  is  determined.  As  was  done  on  the  wing,  spanwise  strips  are 
arranged  on  the  pylon  such  that  angle  0p^_  is  approximately  constant. 

The  strips  are  swept  with  the  pylon  and  have  much  smaller  chord  than 
span  dimension.  In  terms  of  the  coordinate  system  indicated  in  figure  6, 
the  wing  coordinate  system,  the  X-coordinates  of  the  narrow  strip  are 
given  by 


X,  = X + Z,  tan  X 
1 p l 


Pi 


(22) 


where  sweep  angle  Xp^  is  measured  from  the  positive  (downward)  Z-axis, 

positive  for  forward  sweep.  With  the  assumption  of  narrow  strips  and 

for  a pylon  located  at  Yx  = Y , equation  (21)  simplifies  to 

P 


A 0 _ 


tan  0 AX 
P 


-Es-i  f 

27 T J 


dZ 


V(Xp  + tan  *p  - X)2  + (Yd  - Y)2  + (z[  - Z) 

(23) 


Zb  ’ Pc  1 Pi  ' ' P 


where  AX.  ~ X - X . Since  tan  0D,  = tan  (180°  - 6 ) ~ -0  and 
l a d ^t  p p 

Xp^  = -ip  ^ , see  figure  6,  the  velocity  potential  due  to  one  strip  on  the 
pylon  is  then  given  by 


A£  = 


2 

00 


0AXi  r r— 

— 2 — cos  lf,n  — . / (X 


- X - Z 


tan  f.)2  + (Y  - Y)2  + (Z  - Z): 
i p a 


+ — - (X  - X)  sin  ib . - Z cos  ib . 

cos  0^  *1 


- in 


"\/(XPc  " X " Zb  tan  ^i)2  + (Yp  - Y)2  + (Zb  - Z) 


+ rrr~ 7 (X  - X)  sin  f.-Z  cos  f. 

cos  ill-  p Y 1 

1 C 


(24) 
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Differentiating  with  respect  to  X,Y,Z  yields  the  perturbation  velocities 
induced  by  one  spanwise  source  strip  on  the  pylon.  The  results  are 


a ^AX. 

= -ir  cos  ,Ai 

oo 


"N 


"*-X  4*  X + Z tan  ip . + F sin  ip . -X  + X + Z,  tan  ip . 4-  H sin  ip  T 

p a p b *1  yi 


EF 


GH 


a 9 AX  . 

Av  i 


f-Y  4*  Y -Y  + Y 
V„  = 27T  COS  \ 


Aw  0AXi 
— = *”2”  COS  *i  ' 

OO 


where 


-Z  4-  Z - F cos  ib.  -Z,  4-  z - H cos  ib. 

a r l b r l 


EF 


GH 


(25) 


= -\/(X  - X - Z tan  ^.)2  + (Y  - Y) 2 + (Z  - 

V Pc  a 1 p a 


"\ 


Z)2  + 


cos  Ip  ^ 


- (X  - X)  sin  ip . - Z cos  ip . 


F « -^T - X - Za  tan  fL)s  + (Yp  - Y)2  + (Z0  - Z)2 


) (26) 
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G = -v /(X_  - X - Zb  tan  + (Y  - Y)2  + (Zb  - Z)2  + 


cos 


- (X  - X)  sin  tp ^ - Z cos 


H = -^oT  - X - Zb  tan  ^)2  + (Yp  - Y)2  + (Zfe  - Z) 
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At  a field  point  (X,Y,Z)  the  perturbation  velocity  components  induced  by 
pylon  thickness  are  then  given  by  equations  (25)  summed  over  all  pylon 
source  strips  on  both  left  and  right  pylons.  The  specification  of  the 
Yp  coordinate  determines  which  pylon  is  accounted  for. 

In  connection  with  the  flow  tangency  boundary  conditions,  equations 
(6)  and  (7),  perturbation  velocities  induced  by  the  thickness  envelopes 
of  the  wing  and  pylons  are  computed  at  the  control  points  associated 
with  the  vortex  lattice  area  elements  on  the  wing-pylon  combination.  In 
general,  these  velocities  are  obtained  by  adding  the  result  of  summing  equa- 
tions (19)  over  all  wing  source  strips  to  the  result  of  summing  equations  (25) 
over  all  pylon  source  strips.  The  final  sums  are  then  added  to  the  set  of 
externally  induced  perturbation  velocities  u^/V^,  v./Va,  w1/v00  that 
appear  in  the  flow  tangency  boundary  conditions.  In  this  way,  all  inter- 
ference effects  induced  by  the  wing  on  the  pylon  and  vice  versa  are 
accounted  for  and  will  be  reflected  in  the  resulting  circulation  strength 
distribution  on  the  wing-pylon  combination. 

4.3  TER  or  MER  Rack  Flow  Model 

The  remaining  aircraft  component  to  be  modeled  is  the  rack  used  for 
multiple  carriage  of  stores  such  as  a triple  ejector  rack,  TER,  or  a 
multiple  ejector  rack,  MER.  A wind-tunnel  model  of  a TER  rack  is  shown 
in  figure  7.  This  can  be  approximated  by  an  axisymmetric  body  to  which 
are  attached  three  short  pylons.  The  model  used  in  the  present  work  for 
modeling  racks  of  this  type  is  one  which  accounts  only  for  the  volume 
distribution  associated  with  the  axisymmetric  body  portion  of  the  rack. 

This  is  represented  by  a distribution  of  three-dimensional  point  sources 
along  the  longitudinal  axis.  The  source  strengths  and  perturbation 
velocities  are  determined  in  the  same  manner  as  described  in 
section  4.1  for  the  fuselage  and  stores. 

The  short  pylons  have  been  neglected  in  the  present  work.  Their 
effect  could  be  included,  however,  it  has  not  been  done.  The  importance 
of  neglecting  these  short  pylons  will  be  discussed  in  section  7.2. 

5.  FORCE  AND  MOMENT  CALCULATION 

The  calculation  of  the  trajectory  of  a store  dropped  from  an  air- 
craft requires  a knowledge  of  the  forces  and  moments  acting  on  the  store 
at  each  point  in  the  trajectory.  In  order  to  calculate  these  forces  and 
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moments,  the  nonuniform  velocity  field  in  which  the  store  is  immersed,  as 
seen  by  the  store,  must  be  determined  at  each  point  in  time.  This  field 
includes  the  free-stream  velocity,  the  perturbation  velocities  induced  by 
the  parent  aircraft,  and  the  angular  velocities  due  to  the  store's  pitch, 
yaw,  and  roll  motions.  The  distribution  along  the  body  axis  is  required 
in  order  to  calculate  the  forces  and  moments  acting  on  the  body  and,  if 
an  empennage  is  present,  the  distribution  over  the  tail  fin  surfaces  is 
required  for  the  empennage  force  and  moment  calculation. 

This  section  will  first  describe  the  calculation  of  the  velocity 
field  utilizing  the  flow  models  of  section  4.  Following  this  the 
methods  employed  to  calculate  the  body  and  empennage  forces  and  moments 
from  this  nonuniform  flow  field  will  be  presented. 

5.1  Calculation  of  the  Velocity  Field  Including  Damping 

The  xg,  ys,  zg  coordinate  system  to  be  used  in  the  force  and 
moment  calculations  is  shown  in  figure  8.  It  is  fixed  in  the  store  with 
the  origin  at  the  store  nose  with  the  velocities  Ug,  Vg,  and 
positive  in  the  directions  shown.  These  are  the  total  velocities  as 
seen  by  a point  on  the  store  and  each  is  composed  of  a free-stream  component, 
a perturbation  component  due  to  the  parent  aircraft,  and  a damping 
component  due  to  the  store's  rotational  motion. 

Before  writing  the  expressions  for  these  velocities  let  us  examine 
figure  9 which  shows  the  coordinate  systems  used  in  the  trajectory 
calculation.  The  equations  of  motion  will  be  presented  in  section  6. 

Two  coordinate  systems  are  shown  in  figure  9.  The  £,r),£  system  is  an 
inertial  system  fixed  in  the  aircraft  which  is  in  rectilinear  flight  at 
uniform  velocity.  The  positive  £ axis  is  forward  along  the  fuselage 
longitudinal  axis,  the  r]  axis  is  positive  laterally  to  the  right,  and 
the  £ axis  is  positive  downward  perpendicular  to  the  fuselage  longi- 
tudinal axis.  The  origin  of  this  system  is  fixed  in  the  fuselage  nose. 

The  x,  y,  z coordinate  system  is  fixed  in  the  store  with  the 
origin  at  the  store  moment  center.  Referring  back  to  figure  8,  the  x- 
axis  lies  along  the  store  longitudinal  axis  but  is  positive  forward. 

The  y-axis  is  parallel  to  the  yg  axis  and  positive  in  the  same 
direction  and  the  z axis  is  parallel  to  the  zg  axis  but  positive 
downward.  The  x,  y,  z system  shown  in  figure  8 is  a right-hand  system 
and  the  store  angular  velocity  components  about  the  three  axes  are  p,  j, 
and  r respectively. 
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In  figure  9 the  origins  of  the  two  coordinate  systems  have  been 
drawn  to  coincide  in  order  to  show  the  angles  used  to  determine  the  orien- 
tation of  the  store  with  respect  to  the  inertial  (£,rj,£)  axes.  The 
system  of  angles  chosen  consists  of  three  rotations  in  the  yaw,  Y; 
pitch,  6;  and  roll,  $>;  sequence  (ref.  15).  The  two  coordinate  systems 
are  then  related  by 
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The  three  velocity  components,  Us,  Vg,  and  Wg,  are 
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The  first  term  on  the  right  side  of  the  equal  sign  in  each  of  the  above 
expressions  is  the  component  of  the  free-stream  velocity  vector.  If  the 
parent  aircraft  is  moving  at  velocity,  V , and  flying  at  angle  of  attack, 
ct^,  then  the  velocity  of  the  inertial  coordinate  system  relative  to  a 
point  in  space  is 
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(30) 


Vco  = voo  cos  afef  + Vco  sin  afer 

The  velocity  of  the  store  moment  center  relative  to  the  moving  inertial 
system  is 
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so  that  the  velocity  of  the  store  relative  to  a point  in  space  is 


V = \f  + V 
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Since  the  velocity  components  needed  in  equation  (29)  are  those  as  seen 
by  a point  fixed  in  the  store  and  moving  with  the  store,  equation  (32) 
becomes 
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The  components  in  the  x,  y,  z coordinate  system  of  figure  8 are  then 
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where  is  the  transpose,  or  inverse  since  £aJ  is  orthogonal,  of 

the  matrix  given  in  equation  (28) . These  velocity  components  are  related 
to  those  in  the  coordinate  system  of  figure  8 by 
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The  second  term  in  the  three  expressions  of  equation  (29)  is  the 
perturbation  velocity  term.  This  is  comprised  of  the  perturbation  velo- 
cities induced  by  the  fuselage,  wing,  pylon,  rack  and  other  stores  and 
is  calculated  by  the  methods  presented  in  section  4.  It  is  to  be  remem- 
bered that  these  methods  apply  to  the  equivalent  incompressible  configu- 
ration. At  each  point  in  the  trajectory  the  points  at  which  the  velo- 
cities are  required  must  be  located  in  the  incompressible  space.  The 
perturbation  velocities  induced  by  the  fuselage,  rack,  and  store  volumes 
are  calculated  using  equations  (4)  and  (5) , those  induced  by  the  wing- 
pylon  vortex  lattice  are  determined  by  equation  (8)  by  summing  over 
all  vortices,  and  those  induced  by  the  wing  and  pylon  thickness  distribution 
are  obtained  from  equations  (19)  and  (25)  by  summing  over  all  thickness 
strips.  The  velocities  calculated  using  these  equations  are  in  the  coor- 
dinate system  of  each  aircraft  component.  Prior  to  transforming  these 
back  to  the  compressible  space  they  must  be  summed  up  in  the  fuselage 
coordinate  system  since  this  is  the  coordinate  system  in  which  the 
compressibility  correction  was  applied.  This  system  is  the  (£,r],0 
system  of  figure  9. 

Let  u^  , v' , and  w£  be  the  sums  in  the  £,11,  and  £ directions 
respectively.  Then  from  equation  (2)  the  compressible  perturbation 
velocities  are 


(36) 
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The  velocity  components  in  the  x,  y,  and  z directions  of  figure  8 
are 
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and,  finally,  the  components  in  the  coordinate  system  of  figure  8 are 


u = -u,  v = v,  w = -w 
s 9 s 9 s 


(38) 


In  calculating  the  velocity  field  for  the  case  of  a store  released 
from  under  the  fuselage  the  presence  of  the  wing  is  ignored  for  the 
reasons  discussed  in  section  5.1  of  reference  1.  The  arguments  presented 
there  apply  to  the  present  case  even  though  the  configuration  is  more 
complicated.  Fuselage  angle  of  attack  effects  are  included  for  this 
case  and  given  by 


w ' = V sin  a, 1 

^ 00  f 


(39) 


where  wK ' is  the  velocity  normal  to  the  fuselage  axis  in  the  incom- 
pressible space  a distance  from  the  axis,  a is  the  local  fuselage 

radius,  and  a^*  is  the  fuselage  angle  of  attack  in  the  incompressible 
space.  This  velocity  field  is  added  to  the  pylon-induced  and  volume- 
induced  velocity  fields.  Equation  (39)  is  obtained  from  the  complex 
potential  for  a circular  cylinder  in  uniform  flow  given  on  page  29  of 
reference  16.  For  stores  under  the  wing  this  velocity  is  not  included 
in  the  perturbation  field  for  reasons  which  are  also  discussed  in 
section  5.1  of  reference  1. 


The  last  terms  on  the  right  side  of  the  expressions  for  Vg  and 
in  equation  (29)  are  the  yaw  and  pitch  damping  terms.  The  positive 
directions  of  the  rotational  velocities  are  shown  in  figure  9.  it  can 
be  seen  that  r is  positive  for  nose-to-the-right  motion  and  q is 
positive  for  nose-up  motion.  The  means  of  accounting  for  roll  damping 
will  be  discussed  in  the  section  describing  the  empennage  force  and 
moment  calculation,  section  5.3.  The  rotational  velocities,  p,  q,  and  r, 
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are  determined  as  a function  of  time  during  the  integration  of  the  equations 
of  motion. 

The  velocities  to  be  used  in  the  force  and  moment  calculations  are 
made  dimensionless  by  the  store  free-stream  velocity.  Thus,  equation  (29) 
becomes 
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where  from  equation  (32)  or  (33) 

Voo  = RVco  COS  af  + £)2  + + (Voo  sin  af  + C)?]1//2  (41) 

s L 

At  any  point  in  time  during  the  trajectory  the  forces  and  moments  are 
determined  by  removing  the  store  from  the  flow  field,  determining  the 
velocities  at  a series  of  points  which  the  store  longitudinal  axis  and 
tail  fins  occupied,  and  then  immersing  the  store  in  this  flow  field. 

5.2  Body  Forces  and  Moments 

The  store  body  force  and  moment  calculations  have  not  changed 
significantly  since  reference  1.  Slender-body  theory  is  the  basis  for 
calculating  the  body  normal-force  and  side-force  distributions  which 
yield  forces  and  moments.  A simplified  method  of  calculating  the  buoyant 
force  has  been  introduced  which  eliminates  the  necessity  of  calculating 
the  lateral  and  vertical  pressure  gradients  and  speeds  the  computation. 

As  in  reference  1,  when  large  values  of  the  combined  angle  resulting  from 
upwash  and  sidewash  occur,  such  that  the  body  boundary  layer  separates, 
the  slender-body  theory  calculation  is  not  continued.  Simple  viscous 
crossflow  theory  is  used  downstream  of  the  separation  location. 
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5.2.1  Buoyancy  forces  and  moments 


In  the  work  of  reference  1 the  buoyancy  normal-force  distribution 
was  determined  from  the  vertical  pressure  gradient  acting  on  the  body. 

This  required  the  calculation  of  the  velocity  field  not  only  along  the 
body  longitudinal  axis  but  also  along  the  upper  and  lower  surfaces  of 
the  body.  In  the  present  work  the  buoyancy  side-force  distribution  is 
also  required  which  would  necessitate  calculating  the  velocities  along 
the  two  sides  of  the  store  as  well.  Since  calculating  the  velocity  field 
at  five  points  at  each  axial  station  would  be  very  time  consuming, 
a simplified  method  for  calculating  the  buoyancy  forces  and  moments  has 
been  adopted  which  yields  comparable  accuracy. 

This  simplified  method  was  discussed  in  reference  1.  There  it  is 
shown  that  if  it  is  assumed  that  the  flow  in  planes  perpendicular  to  the 
axis  of  the  store  obeys  Laplace's  equation  and  that  the  upwash  and 
sidewash  velocities  vary  along  the  store  length,  and  the  potential  is 
constructed  on  this  basis,  then  integration  of  the  body  pressures  obtained 
from  the  unsteady  Bernoulli  equation  yield  both  the  buoyancy  forces  and 
the  slender-body  forces.  The  final  expression  for  the  buoyancy  normal- 
force  coefficient  is  given  by  equation  (63)  of  reference  1.  The  derivation 
is  carried  out  in  Appendix  I of  that  reference  with  the  final  expression 
which  contains  both  the  buoyancy  and  slender-body  terms  given  by 
equation  (1-14) . 

The  expressions  which  have  been  used  in  the  present  work  to  calcu- 
late the  buoyancy  normal-force  and  side-force  coefficients,  derived  in 
this  manner,  are 
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The  expressions  for  the  pitching-moment  and  yawing-moment  coefficients, 
with  the  moment  taken  about  x (see  fig.  8) , are 
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5.2.2  S lender-body  forces  and  moments 

The  slender-body  forces  and  moments  are  determined  in  the  same 
manner  as  described  in  section  6.2.1  of  reference  1.  The  derivation  will 
not  be  repeated  here.  The  expressions  for  the  forces  and  moments,  from 
that  reference  are 
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The  upper  limit  on  the  above  integrals  is  the  assumed  separation  location, 
that  is,  the  point  at  which  viscous  forces  become  important. 

No  definite  rules  for  the  selection  of  a value  of  x can  be 

s , o 

given.  For  stores  with  cylindrical  afterbodies  xg  o should  probably 
be  taken  as  l,  , the  store  base.  For  stores  with  boattailed  afterbodies 
xs  Q should  probably  be  less  than  {.  . A further  discussion  is  given 
in  sections  6.2  and  8.2  of  reference  1.  In  that  discussion  the  work  of 
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Hopkins  (ref.  17)  is  mentioned.  In  this  work  he  correlated,  for  bodies 
of  revolution  in  uniform  flow,  the  maximum  extent  of  applicability  of 
potential  flow  as  a function  of  the  location  of  the  maximum  negative  rate 
of  change  of  body  cross-sectional  area.  From  that  reference 


xs  o xs  1 

= 0.378  + 0.527 


(50) 


where  xg  Q/ £ s is  the  location  on  the  body  where  potential  flow  is  no 
longer  applicable  and  viscous  forces  are  important  and  xg  £ is  the 
location  on  the  body  of  the  maximum  negative  rate  of  change  of  cross- 
sectional  area.  Whether  or  not  this  correlation  applies  to  nonuniform 
flow  is  not  known  but  it  can  be  used  in  estimating  the  value  of  xg  Q. 

5.2.3  Viscous  crossflow  forces  and  moments 

From  the  assumed  separation  location  to  the  base  of  the  store  a 
viscous  crossflow  calculation  is  used  in  place  of  the  slender  body  calcu- 
lation. The  expressions  for  normal  force  and  pitching  moment  are  derived 
in  section  6.2.2  of  reference  1.  The  expressions  for  side  force  and 
yawing  moment  are  derived  in  a similar  manner.  The  expressions  for  the 
two  forces  and  two  moments  are 
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where 
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and  0^^  is  the  section-drag  coefficient  of  a cylinder  normal  to  the 
free  stream.  That  is 
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5.3  Empennage  Forces  and  Moments 

A method  of  calculating  the  forces  and  moments  associated  with 
planar  or  cruciform  empennages  will  now  be  presented.  It  has  been  noted 
from  experimental  data  that  the  variation  in  flow  velocity  normal  to  the 
empennage  panels  varies  very  little  in  the  chordwise  direction,  for  fins  with 
small  chords,  as  compared  to  the  spanwise  direction.  Accordingly,  the  effec- 
tive shape  of  the  empennage  panels  is  essentially  that  of  pure  twist  with  zero 
camber.  For  such  panels,  it  is  possible  to  calculate  the  empennage  forces  using 
reverse-flow  theorems.  The  net  force  generated  by  a panel,  including  that 
carried  over  onto  the  store  body,  can  be  obtained  in  terms  of  an  integral, 
across  the  span  of  the  panel,  of  the  product  of  the  twist  angle  times 
an  influence  function  which  corresponds  to  span  loading  of  an  associated 
empennage  in  uniform  reverse  flow.  Influence  functions  have  been  derived 
for  calculating  empennage  normal  force,  side  force,  pitching  moment, 
yawing  moment,  and  rolling  moment  for  planar  and  cruciform  empennages. 

For  cruciform  empennages  panel-panel  interference  has  been  accounted  for. 
Complete  mathematical  details  of  the  method  are  contained  in  Appendix  I. 

A crossflow  plane  in  the  empennage  region  is  shown  in  figure  10 
for  a cruciform  empennage.  The  tail  fins  are  numbered  Fx , F2,  F3,  and  F4 
such  that,  for  a planar  empennage,  fins  F3  and  F4  are  not  present.  This 
crossflow  plane  is  the  plane  in  which  the  empennage  forces  are  assumed  to 
act  and  will  be  taken  as  that  which  contains  the  quarter  chord  of  the 
mean  aerodynamic  chord  of  the  exposed  panels.  In  the  present  work  all  of 
the  fins  will  be  assumed  to  have  the  same  planform. 

The  x^,  y^,  z^  x coordinate  system  is  fixed  in  the  tail  fins  such 
that  the  x^,  y^  plane  contains  the  horizontal  fins  and,  in  the  case 
of  a cruciform  empennage,  the  xf,  zf  plane  contains  the  vertical  fins. 
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This  coordinate  system  is  rolled  an  angle  0^  from  the  xg,  yg,  zg 
store  body  coordinate  system  of  figure  8.  Both  coordinate  systems  move 
with  the  store  during  the  trajectory.  In  the  trajectory  calculation  the 
xg,  zg  plane  is  assumed  to  be  parallel  to  the  plane  in  the  inertial 

system  (see  fig.  9)  at  time  t = 0 seconds.  Thus,  the  specification  of 
the  angle  0^  allows  a trajectory  to  be  started  with  the  fins  in  other 
than  the  vertical  and  horizontal  positions. 

Also  shown  in  figure  10  are  the  velocities  required  for  the  empennage 
force  and  moment  calculation.  The  arrows  indicate  the  positive  directions. 
The  velocities  on  each  fin  are  required  at  n equally  spaced  points 
between  the  body  radius  and  the  tip  of  the  panel.  The  coordinates  of 
the  points  are  (i  = l,...n) 
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The  x^  location  of  all  points  is  the  same,  that  is,  the  axial  location 

of  the  quarter  chord  of  the  mean  aerodynamic  chord,  in  figure  8.  In 

calculating  the  velocities  normal  to  the  fins  the  coordinates  of  the  points 

given  by  equation  (56)  are  first  transformed  to  the  store  body  coordinate 

system,  the  xg,  yg,  zg  system,  and  then  the  velocities  in  this  system 

are  calculated  as  described  in  section  5.1.  These  velocities  are  given 

by  equation  (29)  and  include  the  store  translational  motion  as  well  as 

the  yaw  and  pitch  rotations.  The  velocities  normal  to  the  fins  are  then 

obtained  by  summing  up  the  components  of  W and  V normal  to  the  fins. 

s s 
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W,  = W cos  <b£  + V sin  <f>  - 

1 s s 

W^>  = W cos  (hr  + V sin  0- 

2 s s 

V1  = -W  Sin  <br-  + V COS  0- 

1 s s 

Vo  = -W  sin  0-  + V cos  0-^ 

2 s f s f 


(57) 


The  equations  presented  in  Appendix  I determine  the  forces  and 
moments  produced  by  the  empennage  in  the  x^,  y^,  coordinate  system 

of  figure  10.  The  force  given  by  equation  (1-13)  is  positive  in  the  z^ 
direction  and  the  force  given  by  equation  (1-18)  is  positive  in  the  y^ 
direction.  The  moment  given  by  equation  (1-21)  is  positive  nose  up  about 
the  y^  axis  and  that  of  equation  (1-22)  is  positive  nose  to  the  right 
about  the  z^  axis.  Since  the  x^  and  xg  axes  coincide  the  rolling 
moment  expressions,  equations  (1-30)  and  (1-52) , apply  to  both  coordinate 
systems  and  the  positive  sense  is  right  wing  down. 

In  the  store  body  coordinate  system,  which  is  the  system  the  forces 
and  moments  must  be  in  for  the  trajectory  calculation,  the  normal-force 
and  side-force  coefficients  are,  from  equations  (1-13)  and  (1-18) , 
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Similarly,  from  equations  (1-21)  and  (1-22) , the  pitching-moment  and 
yawing-moment  coefficients  are 
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The  rolling-moment  coefficient  for  a planar  empennage  is  given  by  equa- 
tion (1-30) 


(CJ  - (CJ 

1 E 1 H 


(62) 


and  that  for  a cruciform  empennage  by  equation  (1-52) 
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HV 
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6.  CALCULATION  OF  STORE  TRAJECTORIES 

In  the  preceding  sections  of  this  report,  methods  have  been  presented 
for  modeling  the  various  aircraft  components,  for  calculating  the  non- 
uniform  flow  field  due  to  the  aircraft,  and  for  calculating  the  forces 
and  moments  acting  on  a store  immersed  in  this  flow  field  at  any  location 
relative  to  the  aircraft.  Utilizing  these  tools  a computer  program  has 
been  written  which  integrates  the  equations  of  motion  of  the  store  to 
determine  its  location  and  orientation  relative  to  the  parent  aircraft  as 
a function  of  time.  In  the  work  of  reference  1,  this  motion  was  restricted 
to  three  degrees  of  freedom.  In  the  present  work  the  equations  of  motion 
allow  six  degrees  of  freedom. 

In  this  section  of  the  report  the  equations  of  motion  will  first  be 
discussed  and  then  the  computer  program  will  be  described. 

6.1  Equations  of  Motion 

The  complete  derivation  of  the  equations  of  motion  of  a rigid  body 
with  mass  and  inertia  asymmetries  is  presented  in  Appendix  II.  Generally 
the  body,  a store  in  the  present  application,  will  have  axes  of  geometri- 
cal symmetry  about  which  the  store  forces  and  moments  are  determined.  By 
mass  asymmetry  it  is  meant  that  the  store  center  of  mass  does  not  have  to 
lie  at  the  origin  of  this  coordinate  system.  By  inertia  asymmetry  it  is 
meant  that  the  principal  axes  of  inertia  of  the  store  do  not  coincide 


with  the  geometric  axes  of  symmetry  so  that  the  tensor  of  inertia  possesses 
off-diagonal  terms,  the  products  of  inertia. 

In  the  derivation  of  Appendix  II  the  two  coordinate  systems  discussed 
earlier  in  conjunction  with  the  velocity  field  calculation,  section  5.1, 
are  used.  They  are  shown  in  figure  9.  The  x,  y,  z coordinate  system 
is  fixed  in  the  store  and  rotates  with  the  store.  The  x-axis  lies  along 
the  store  longitudinal  axis  and  is  positive  forward,  the  y-axis  is  posi- 
tive to  the  right,  and  the  z-axis  is  positive  downward.  These  axes  coin- 
cide with  the  geometric  axes  of  symmetry  with  the  origin  at  the  point  about 
which  the  moments  and  products  of  inertia  and  the  aerodynamic  moments  are 
calculated.  The  equations  of  motion  are  derived  in  this  coordinate  system 
since,  by  allowing  the  coordinate  system  to  rotate  with  the  store,  the 
time  derivatives  of  the  moments  and  products  of  inertia  do  not  appear  in 
the  rotational  equations  of  motion.  The  equations  of  motion  in  this 
coordinate  system  are  given  by  equations  (11-15)  and  (11-40)  of 
Appendix  II. 

Since  the  store  coordinate  system  in  which  these  equations  of  motion 
are  written  does  not  allow  the  store  position  and  orientation  to  be  deter- 
mined relative  to  the  parent  aircraft,  a coordinate  system  is  introduced 
which  is  fixed  in  the  fuselage  nose,  the  £ , rj , £ system  of  figure  9.  In 
the  present  work  it  is  assumed  that  the  parent  aircraft  is  flying  at 
constant  velocity,  constant  angle  of  attack,  and  constant  flight  path 
angle  relative  to  the  horizontal  so  that  this  system  is  a non-rotating 
system.  As  such,  it  can  be  regarded  as  an  inertial  system  and  the  store 
motion  relative  to  this  moving  system  calculated. 

The  £,  rj,  £ system  if  fixed  in  the  fuselage  nose  with  £ forward 
along  the  fuselage  longitudinal  axis,  rj  laterally  to  the  right,  and  £ 
vertically  downward.  In  order  to  determine  the  orientation  of  the  store 
with  respect  to  these  axes,  a system  of  angles  must  be  introduced.  These 
angles  are  also  shown  in  figure  9 and  consist  of  three  rotations  in  the 
yaw,  pitch,  9;  roll,  <J>;  sequence.  In  order  to  determine  the  time 
histories  of  these  angles  differential  equations  expressing  them  as  a 
function  of  the  store  rotational  velocities,  p,  q,  and  r,  are  required. 
These  are  given  by  equation  (II-l)  of  Appendix  II.  The  rotational 
velocities,  p,  q,  and  r,  are  shown  in  figure  9.  The  velocity  about  the 
x-axis  is  p and  is  positive  right-wing  down.  Those  about  the  y-  and 
z-axes  are  q and  r respectively  and  are  positive  nose  up  and  nose  to  the 
right. 
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The  final  differential  equations  which  are  to  be  integrated  to 
determine  the  store  position  and  orientation  relative  to  the  inertial 
system  fixed  in  the  fuselage  are  equations  (II-l) , (11-16)  through  (11-18) , 

and  (11-41)  through  (11-43) . The  latter  six  equations  are  obtained  from 
equations  (11-15)  and  (11-40)  by  using  the  direction  cosines  relating 
the  two  coordinate  systems,  equation  (II-2) . In  the  equations,  £,  q,  and 
£ are  the  accelerations  of  the  store  moment  center,  the  point  about 
which  the  moments  and  products  of  inertia  are  taken,  in  the  inertial 
coordinate  system.  The  velocities  in  this  system  are  I*,  r),  and  £ and 
the  quantities  xq,  yQ,  and  zq  are  these  velocities  resolved  into  the 
store  body  coordinate  system.  That  is 


Z 

oJ 


(64) 


where  H is  the  transpose  of  the  matrix  given  by  equation  (II-2)  . 

The  location  of  the  store  center  of  mass  relative  to  the  store  moment 
center  in  store  body  coordinates  is  x,  y,  z.  The  three  rotational  velo- 
cities in  the  store  body  coordinate  system  are  p,  q,  and  r and  the 
corresponding  accelerations  are  p,  q,  and  r.  The  quantities  1^,  I , 

and  I are  the  moments  of  inertia  and  1,1,  and  I are  the 
zz  yz*  xz*  xy 

products  of  inertia.  They  are  defined  by  equation  (11-36) . 

The  remaining  parameters  in  the  equations  of  motion  are  the  store 

mass,  m,  and  the  forces,  F , F , and  F , and  moments,  M , M , and  M , 

x y z x y z 

acting  on  the  store  which  are  positive  in  the  positive  x,  y,  z,  p,  q, 
and  r directions  respectively.  The  three  forces  are 


Fx  = m9x  - ^ SrS0 
s 

Fy  ' ”9y  + VrCy  > 

Fz  - mg2  - <^  SC 
s 


(65) 
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In  these  equations  g , g , and  g are  the  component  of  the  gravitational 

x y z 

acceleration  vector  "g  in  the  store  body  coordinate  system.  Consider  the 
following  sketch 


which  shows  the  fuselage  flying  at  velocity  V^,  angle  of  attack  af, 
and  flight  path  angle  7 The  gravitational  force  acts  normal  to  the 
local  horizontal  so  that  in  the  inertial  coordinate  system 


Then 


g = -g  sin  (af  + ) ef  + g cos  (af  + 7f ) e^ 


V 

g sin  (af  + yf) 

gy 

-or 

0 

1 

N 

* j 

g cos  (af  + Yf) 

(66) 
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The  second  terms  in  equation  (65)  are  the  aerodynamic  forces  acting 
on  the  store.  The  axial-force  coefficient,  C^,  is  specified  and  is  posi- 
tive in  the  negative  x direction.  The  sideforce  and  normal  force 
coefficients,  CY  and  CN,  are  calculated  by  the  methods  described  in 
section  5 and  are  positive  in  the  positive  y and  negative  z directions 
respectively.  These  coefficients  are  the  sums  of  all  of  the  components 
calculated  in  section  5.  For  example 
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The  parameters  a and  SD  are  the  dynamic  pressure  and  reference  area 
■^s  K 

used  in  nondimens ionalizing  the  forces.  In  the  present  work  they  are 
taken  as  the  store  free-stream  dynamic  pressure  and  the  store  maximum 
cros  s-sectional  area 


^00  2 P 00^00 

s s 


(69) 
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The  density  pm  is  assumed  constant  at  the  value  for  the  aircraft  flight 
altitude  at  time  t = 0 and  the  store  free-stream  velocity  is  given  by 
equation  (41) . 

The  three  moments  about  the  store  moment  center  are 


— 


MX  = Vr0*  + m(gzy  - gyz) 


M = a Sni_C  + m (g  z - g x) 
y R R m yx  ^z  ' / 

J C 


(71) 


M 


= a S^IC  + m(g  x - g y) 
z ^ R R n v^y  V J 


The  first  term  on  the  right-hand  side  of  each  of  the  above  expressions  is 

the  aerodynamic  moment.  The  moment  coefficients,  C^,  C^,  and  Cn,  are 

calculated  by  the  methods  of  section  5.  The  reference  length,  l , is 

R 

taken  as  the  maximum  store  diameter 


^R  ^amax 


(72) 


The  second  terms  are  moments  produced  by  the  gravitational  force  when 
the  location  of  the  store  center  of  mass  does  not  coincide  with  the 
center  of  moments. 

In  order  to  integrate  the  equations  of  motion  given  by  equations 
(11-16)  through  (11-18)  and  (11-41)  through  (11-43)  and  the  equations 
determining  the  angular  orientation,  equation  (II-l)  , initial  conditions 
of  the  variables  must  be  specified.  These  are  the  position  of  the 
store  center  of  moments  (£,  t},  £)  , the  translational  velocity  of  this 
point  (£,  ij , C)  > the  angular  velocities  about  the  three  store  axes 
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p,  q,  r)  , and  the  initial  orientation  angles  (¥,  0,  $)  . With  these 
specified  the  equations  can  be  integrated  with  the  aerodynamic  forces  and 
moments  being  recalculated  at  each  point  in  the  trajectory. 

6.2  Description  of  Computer  Program 

The  computer  program  which  has  been  written  to  calculate  the  trajectory 
of  a store  ejected  from  an  aircraft  consists  of  two  programs,  an  axi- 
symmetric  source  distribution  program  and  a trajectory  program.  The  calcu- 
lation of  the  source  distributions  to  represent  the  various  axisymmetric 
bodies  has  been  kept  separate  for  two  reasons.  First,  a number  of  runs 
may  have  to  be  made  with  the  program  to  obtain  a source  distribution  which 
adequately  reproduces  the  shape  of  the  body.  The  second  reason  is  that 
once  the  source  distribution  is  obtained  for  a specific  Mach  number  it 
need  not  be  recalculated  since  for  a given  body  shape  the  distribution  is 
only  a function  of  Mach  number. 

For  a given  aircraft-stores  combination  and  Mach  number  source 
distributions  must  be  obtained  for  the  fuselage,  each  different  store 
shape,  and  the  ejection  rack  if  one  is  present.  The  various  shapes  are 
specified  by  a series  of  segmented  polynomials  and  the  program  calculates 
and  prints  the  source  locations  and  strengths.  The  locations  are  deter- 
mined by  the  program  from  two  input  parameters,  the  location  of  the  source 
nearest  the  nose  and  the  distance  between  sources  which  is  specified  as 
a fraction  of  the  local  body  radius. 

The  trajectory  program  consists  of  three  main  sections.  The  first 
section  reads  in  the  input  data  and  calculates  quantities  from  these  data, 
the  second  section  solves  for  the  vorticity  distribution  which  represents 
the  wing-pylon  loading,  and  the  third  section  calculates  the  trajectory. 

The  input  data  section  of  the  program  reads  in  the  following  infor- 
mation: 

(1)  Aircraft  flight  conditions 

(2)  Indices  specifying  what  aircraft  components  are  present 

(3)  Fuselage  data 

(4)  Wing  data 

(5)  Pylon  data 

(6)  Rack  data 

(7)  Store  data 
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The  aircraft  flight  conditions  are  the  angle  of  attack,  flight  path 
angle,  Mach  number,  free-stream  air  density,  and  flight  velocity. 

Item  (2)  consists  of  four  indices.  The  first  three  specify  whether 
or  not  the  fuselage,  pylon,  and  rack  are  present  and  the  fourth  specifies 
the  number  of  stores.  The  program  in  its  present  form  will  handle  a 
maximum  of  ten  stores. 

The  fuselage  input  data  consist  of  the  length,  maximum  radius,  and 
source  distribution. 

The  wing  input  data  locate  the  wing  relative  to  the  fuselage,  inclu- 
ding its  incidence,  and  supply  information  required  to  lay  out  the  vorticity 
and  thickness  distributions.  The  data  also  include  the  twist  and  camber 
distribution  and  the  slope  distribution  of  the  thickness  envelope. 

The  pylon  input  data  locate  the  pylon,  from  which  the  store  is  to 
be  ejected,  laterally  relative  to  the  fuselage  longitudinal  axis  and 
longitudinally  relative  to  the  leading  edge  of  the  local  wing  chord. 

As  for  the  wing,  data  required  to  lay  out  the  vorticity  and  thickness 
distributions  are  input  together  with  the  slope  distribution  of  the 
thickness  envelope.  The  pylon  is  assumed  to  have  zero  twist  and  camber. 

The  rack  is  assumed  to  be  immediately  below  the  pylon.  The  input 
data  for  the  rack  locate  its  nose  relative  to  the  local  wing  chord  and 
specify  its  incidence  relative  to  the  wing  root  chord.  In  addition  the 
rack  length  and  maximum  radius  are  input. 

For  each  store  present  on  the  aircraft,  data  are  input  which  assign 
it  a store  number  and  specify  a shape  number,  the  length,  and  the  maximum 
radius.  Each  store  is  located  by  specifying  its  lateral  position  relative 
to  the  wing  root  chord  and  the  longitudinal  and  vertical  position  of  the 
store  nose  relative  to  the  wing  chord  immediately  above  the  store.  The 
incidence  of  the  store  relative  to  the  wing  root  chord  is  also  specified. 

A source  distribution  is  input  for  each  different  store  shape. 

The  second  main  part  of  the  program  calculates  the  wing-pylon  vor- 
ticity distribution.  This  is  done  using  equations  (6)  and  (7)  of 
section  4.2.2.  The  coefficient  matrix  is  first  calculated  and  then  the 
right-hand  side  is  determined.  The  strengths  of  the  vortices  are  then 
calculated  by  solving  the  set  of  simultaneous  equations. 
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The  third  and  last  main  part  of  the  trajectory  program  is  the 
trajectory  calculation.  This  section  of  the  program  consists  of  the 
following  steps: 

(1)  Input  additional  information  to  describe  ejected  store. 

(2)  Input  empennage  data  if  store  has  one. 

(3)  Initialize  for  trajectory  calculation. 

(4)  Calculate  aerodynamic  forces  and  moments. 

(5)  Calculate  accelerations  and  rates  of  change  of  orientation 

angles. 

(6)  Integrate  equation  of  motion. 

The  additional  data  describing  the  store  to  be  ejected  include  indices 
specifying  which  store  is  to  be  ejected,  the  number  of  segments  the  body 
is  to  be  broken  into  for  the  force  calculation,  the  separation  location, 
and  whether  the  store  has  an  empennage.  Also,  the  store  mass  and  inertia 
characteristics  are  read  in  along  with  the  location  Of  the  point  about 
which  the  aerodynamic  moments  are  to  be  calculated.  This  is  the  point 
about  which  the  moments  and  products  of  inertia  were  taken.  The  location 
of  the  store  center  of  mass  relative  to  this  point  is  also  specified  as 
are  the  store  axial  force  coefficient  and  the  value  of  the  crossflow 
drag  coefficient  to  be  used  in  the  viscous  crossflow  force  and  moment 
calculation. 

Three  other  indices  are  input  which  pertain  to  options  included  in 
the  computer  program.  Provision  has  been  made  to  include  or  exclude  the 
damping  terms  in  the  velocity  field  calculation.  Also,  for  stores  with 
empennages,  rolling  moment  may  or  may  not  be  included  in  the  accleration 
determination.  The  third  option  pertains  to  the  calculation  of  free- 
flight  trajectories  as  opposed  to  captive-store  trajectories  as  obtained 
in  the  wind  tunnel.  In  wind-tunnel  captive-store  testing  it  is  customary 
to  change  the  store  pitch  and  yaw  angles  to  account  for  translational 
motion  only  while  measuring  the  aerodynamic  forces  and  moments.  This 
changes  its  position  in  the  nonuniform  flow  field.  Provision  has  been 
made  in  the  computer  program  to  simulate  this. 

For  a store  with  an  empennage,  additional  quantities  must  be 
specified.  These  data  are  an  index  indicating  whether  the  empennage  is 
planar  or  cruciform,  the  tail-fin  semispan,  the  average  body  radius  in 
the  tail-fin  region,  the  initial  roll  orientation  of  the  fins,  and  the 
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lift-curve  slope  of  the  fins  alone.  In  addition  the  axial  position  at 
which  the  forces  are  assumed  to  act  must  be  specified. 

In  the  trajectory  initialization  certain  store  ejection  conditions 
are  specified.  These  are  an  initial  velocity  in  the  vertical  plane 
perpendicular  to  the  store  longitudinal  axis  and  an  initial  pitching 
velocity.  The  store  lateral  velocity  and  yawing  and  rolling  velocities 
are  initially  zero. 

The  integration  of  the  equations  of  motion  is  done  by  a standard 
numerical  integration  technique  with  the  aerodynamic  forces  and  moments 
calculated  at  each  point  required  by  the  integration  scheme.  The  calcu- 
lation of  the  nonuniform  velocity  field  and  the  resulting  forces  and 
moments  was  described  in  section  5. 

7.  COMPARISONS  WITH  EXPERIMENTAL  DATA 

This  section  of  the  report  will  present  comparisons  with  experi- 
mental data  in  order  to  assess  the  accuracy  of  the  theoretical  methods 
used  to  predict  the  nonuniform  flow  field,  store  body  loading  distribution 
and  forces  and  moments,  empennage  forces,  and  store  trajectories.  In 
addition  MER  interference  effects  are  examined  as  is  the  effect  of  store 
boattail  on  the  loadings  and  forces.  Many  more  comparisons  have  been 
made  with  experimental  data  than  are  presented  here.  The  degree  of 
agreement  to  be  shown,  however,  is  typical  of  that  obtained  for  other 
conditions . 

All  of  the  data  comparison  to  be  shown  use  wind-tunnel  data  from 
references  2 and  3.  The  wind-tunnel  model  used  in  these  investigations 
is  an  idealized  fighter-bomber  aircraft.  The  wing- fuselage  combination 
is  shown  in  figure  11  and  the  pylon,  MER  rack,  and  stores  used  in  con- 
junction with  the  wing-fuselage  are  shown  in  figure  12.  The  TER  rack 
used  was  shown  in  figure  7.  The  large  store  with  cylinderical  afterbody 
shown  in  figure  12(c)  was  tested  singly  with  and  without  empennage  and 
in  the  TER  configuration  without  empennage.  The  large  store  with  boattail 
afterbody  shown  in  figure  12(d)  was  also  tested  singly  with  empennage 
to  evaluate  boattail  effects.  Pressure  instrumented  versions  of  both 
large  stores  without  empennages  were  used  for  obtaining  loading  distri- 
butions. The  forces  and  moments  for  the  boattailed  store  without  empennage 
were  obtained  by  integrating  the  loading  distributions.  In  figure  12(e) 
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is  shown  the  small  boattail  store  which  was  tested  in  the  MER  configuration 
under  the  fuselage  with  and  without  the  empennage.  A table  of  ordinates 
for  the  body  of  this  store  is  given  in  Table  I of  reference  1. 


In  the  discussion  which  follows,  the  various  aircraft  components 


are  designated  as 

w 

wing 

B 

fuselage 

P 

pylon 

S 

c 

large  store  with  cylindrical  afterbody  tested 
without  other  stores 

sb 

large  store  with  boattail  afterbody 

^1  > ^2  > ^3 

large  store  with  cylindrical  afterbody  when 
tested  in  TER  configuration.  The  cyclic  order 
is  shown  below  viewed  from  the  rear 

^1 > Sg , S 3 , S 4 ,S5 ,Sg 

small  boattail  stores  tested  in  MER  configu- 
ration. SijSg,  and  S3  are  the  front  three 
with  the  same  cyclic  order  as  the  TER  con- 
figuration and  S4,Ss,S6  are  the  rear  three 
in  the  corresponding  cyclic  order 

T 

TER  rack 

M 

MER  rack 

Whether  the  empennages  are  present  is  taken  from  the  context.  In 
addition,  the  parent  aircraft  and  the  store  are  at  the  same  angle  of 
attack  in  all  cases  to  be  shown. 

7.1  Flow  Field  and  Loading  Predictions  for  Single  Stores 

The  comparisons  to  be  presented  in  this  section  are  all  for  the 
configuration  shown  in  figure  13.  The  large  store  with  cylindrical 
afterbody,  but  without  the  cruciform  fins  shown  in  figure  12(c),  is 
mounted  at  the  one-third  semispan  location  y/s  = -0.333.  The  pylon  is 
not  shown  in  figure  13  but  the  store  is  shown  in  the  position  it  would 
be  in  if  it  were  attached  to  the  pylon. 
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7.1.1  Effect  of  pylon  on  flow  field  and  attached  store  load  distribution 

The  upwash  angle  ws/v00  and  the  sidewash  angle  V /V^  are  shown 
in  figure  14  for  6°  angle  of  attack  and  Mach  number  0.25  at  the  position 
the  centerline  of  the  attached  store  would  occupy  if  it  were  present  at 
the  one-third  semispan  position  of  the  left  wing  panel.  Data  with  and 
without  the  pylon  present  are  compared  with  theory.  The  sidewash  is 
outboard  everywhere  and  there  is  a small  measured  effect  of  the  pylon 
on  the  sidewash  which  is  less  than  the  predicted  effect.  There  is  a large 
upwash  along  the  store  centerline,  and  its  modification  due  to  the  pylon 
is  well  predicted  by  the  theory. 

The  normal-force  distributions  along  the  store  due  to  the  upwash 
fields  just  described  are  shown  in  figure  15(a).  It  is  seen  that  the 
load  distribution  with  the  pylon  present  is  very  irregular,  and  it  is 
predicted  quite  well.  A number  of  stations  along  the  store  axis  are 
necessary  to  describe  such  a load  distribution  if  it  is  to  be  integrated 
to  obtain  the  force  and  moment.  A minimum  number  of  about  20  stations 
is  required.  It  is  noted  from  the  previous  figure  that  the  upwash  is 
positive  everywhere  and  has  an  average  value  corresponding  to  about  3.5°. 
Behind  the  shoulder  of  the  store  all  the  loading  is  due  to  flow  curvature 
(rate  of  change  of  upwash  with  axial  distance)  since  the  store  after- 
body is  cylindrical.  It  is  seen  that  the  net  lift  is  nearly  zero,  and 
the  influence  of  flow  curvature  is  very  important  in  determining  the 
force  and  moment  on  the  store  in  the  present  attached  position.  The 
effect  of  the  pylon  is  pronounced. 

The  data  shown  in  figure  15(b)  exhibit  a small  influence  of  the 
pylon  on  the  side-force  distribution  for  the  same  conditions  shown  in 
figures  14  and  15(a).  The  small  magnitude  of  the  pylon  effect  is  some- 
what under  predicted. 

7.1.2  Effect  of  store  position  and  angle  of  attack  on  store  forces 

As  the  store  drops  below  its  attached  position  for  the  configuration 
of  figure  13,  the  influence  of  the  pylon  on  the  loading  distributions 
diminishes  rapidly  but  the  gross  forces  do  not  change  greatly.  Consider 
figure  16  which  shows  the  variation  of  store  normal  force  and  side  force 
with  vertical  location  at  a fixed  angle  of  attack  of  6°.  The  distance 
Az  is  zero  for  the  attached  position.  There  is  not  much  change  in  normal 
force  or  side  force  acting  on  the  store  for  several  diameters  beneath  the 
pylon,  and  the  theory  predicts  this  trend  well. 
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In  figure  17,  the  variations  with  angle  of  attack  of  the  store 
normal-force  coefficient  and  side-force  coefficient  are  shown.  The  store 
is  one  diameter  below  the  attached  position  Az/d  « 1.0.  The  variations 
with  angle  of  attack  are  suprisingly  small  and  are  predicted  well  by  the 
theory. 

7.2  Flow-Field  and  Loading  Predictions  for  TER  Configurations 

In  this  section,  the  comparisons  to  be  presented  are  for  the  TER 
configuration  shown  in  figure  18.  The  large  stores  with  cylindrical 
afterbody  and  no  fins  are  used  in  the  TER  grouping.  As  shown  in  figure  18, 
the  stores  are  in  the  carriage  position  when  the  pylon  of  figure  12(a) 
and  the  TER  of  figure  7 are  present.  They  are  not  shown  in  figure  18 
for  clarity.  As  in  the  single  store  case,  the  stores  are  at  the  one- 
third  semispan  location.  For  all  comparisons  to  be  shown  the  Mach  number 
is  0.25. 

7.2.1  Flow-field  predictions 

The  flow  field  has  been  measured  for  a number  of  positions  beneath 
a TER  configuration  with  the  bottom  store  missing.  The  TER  rack  is  longer 
than  the  pylon  and  has  been  approximated  by  an  equivalent  body  or  revo- 
lution in  the  prediction  method.  The  values  of  V /V  and  W /W  are 

s °°  s °° 

shown  in  figure  19  along  the  position  the  axis  of  store  S1  would 

occupy  for  the  attached  store  location,  Az/d  = 0,  and  for  one  diameter 

below  the  attached  location  Az/d  = 1.0.  In  the  figure  title,  WBPTS2S3 

indicates  the  parent  aircraft  configuration  generating  the  flow  field. 

There  is  considerable  difference  in  the  upwash  between  these  vertical 

locations  which  is  well  predicted  by  the  theory.  In  the  theory  the 

effect  of  the  shoulder  stores  S2  and  S3  has  been  predicted  using  three- 

dimensional  source  terms.  Higher  harmonics  such  as  dipoles  fall  off 

more  rapidly  with  distance  from  the  store  than  source  terms  which  dominate 

the  flow  far  away.  It  is  clear  that  neglecting  these  high  harmonics 

yields  fairly  accurate  results  for  both  V /V  and  W /V  in  the 
u s 00  s 00 

s s 

present  instance. 

The  variations  with  angle  of  attack  of  the  sidewash  and  upwash 
angles  are  shown  in  figure  20  for  a position  of  the  store  one  diameter 
beneath  its  attached  position.  Both  the  sidewash  and  upwash  angles  are 
predicted  well  for  a = 0°  and  a = 6°.  The  calculative  method  thus 
predicts  well  the  effects  of  vertical  location  and  angle  of  attack  of  the 
aircraft  on  the  flow  field  acting  on  the  store. 
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7.2.2  Loading  predictions 

The  normal-force  distribution  on  store  Sx  of  the  TER  configuration 
for  the  attached  store  and  for  a location  one  diameter  below  its  attached 
position  are  shown  in  figure  21(a).  The  loading  on  the  attached  store 
is  somewhat  underpredicted  in  the  region  of  the  pylon.  However,  at  one 
diameter  beneath  the  pylon  the  prediction  is  very  good.  Similar  results 
for  the  side-force  distribution  are  shown  in  figure  21(b)  with  the  same 
general  results  for  comparison  between  theory  and  experiment.  The  loading 
theory  utilizes  predicted  flow  field  quantities  and  does  not  include  the 
special  rack  effect  discussed  subsequently. 

The  effect  of  angle  of  attack  on  the  normal-force  distribution  of 
Sr  of  the  TER  configuration  is  shown  in  figure  22(a)  for  a distance  of 
the  store  one  diameter  beneath  its  attached  position.  Experiment  and 
theory  are  in  good  agreement  for  the  angles  of  attack  of  0°  and  6° 
shown  in  the  figure.  Analogous  results  are  shown  in  figure  22(b)  for 
the  side-force  distribution  with  good  comparison  between  experiment  and 
theory. 

The  effect  of  the  rack  on  the  loading  distribution  on  store  Sx  in 
the  attached  position  is  shown  in  figure  23.  The  rack,  shown  in  figure  7, 
can  be  modeled  approximately  by  a body  of  revolution  with  a pylon  beneath 
it.  The  theory  shown  in  figure  23  for  the  rack  effect  is  based  solely 
on  the  source  distribution  corresponding  to  its  equivalent  body  of  revo- 
lution. As  such  it  causes  the  change  in  theoretical  load  distribution 
given  in  figures  23(a),  but  it  has  no  effect  on  the  side-force  distribution 
shown  in  figure  23(b).  Actually,  the  effect  of  the  rack  on  the  normal- 
force  distribution  exhibited  by  figure  23(a)  is  qualitatively  similar  to 
the  effect  of  the  pylon  on  the  normal-force  distribution  on  a single 
store  as  shown  in  figure  15(a).  Accordingly,  the  discrepancy  between 
theory  and  experiment  shown  in  figure  23(a)  can  probably  be  ascribed  to 
neglect  of  the  rack  pylon  in  the  model.  This  effect  is  of  importance 
only  for  attached  loads,  and  is  negligible  for  the  store  one  diameter 
below  its  attached  position.  It  can  be  included  in  the  present  computer 
program,  but  has  not  been  at  this  time.  The  rack  pylon  can  also  have 
an  effect  on  side-force  distribution  as  shown  in  figure  23(b).  The 
theoretical  curves  for  the  two  cases  are  identical  since  the  source  distri- 
bution representing  the  rack  produces  no  side  load. 
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The  loading  distributions  on  store  S2  of  the  TER  configuration 
after  Sx  has  been  released  are  shown  in  figure  24  for  the  attached 
position  and  one  diameter  beneath  the  attached  position.  For  the  attached 
position  there  are  effects  on  normal  force  and  side  force,  which  are  probably 
associated  with  the  secondary  pylon,  not  accounted  for  by  the  theory. 

At  a distance  of  one  diameter  beneath  its  attached  position,  the  experi- 
ment and  theory  are  in  good  agreement  for  both  normal  force  and  side 
force. 

For  a position  of  S2  one  diameter  below  its  attached  position  the 
loading  distributions  are  shown  for  a = 0°  and  a = 6°  in  figure  25. 

The  normal  force  distribution  is  predicted  quite  well  for  both  angles  of 
attack  in  figure  25(a)  and  the  side-force  distribution  predictions  shown 
in  figure  24(b)  exhibit  similar  agreement. 

7.3  MER  Interference  Effects 

The  interference  effects  previously  described  for  a TER  configuration 
also  exist  in  a MER  configuration,  but  in  addition  there  is  interference 
between  the  front  three  stores  and  the  rear  three  stores.  A MER  configu- 
ration has  been  tested  under  the  fuselage  as  shown  in  figure  26.  The 
stores  used  on  the  MER  rack  of  figure  12(b)  are  those  shown  in  figure  12(e). 
Tests  have  been  conducted  for  stores  with  and  without  cruciform  empennages. 
For  all  comparisons  to  be  shown  the  Mach  number  is  0.40. 

In  figure  27,  the  squares  indicate  the  influence  of  the  front  three 
stores,  SlyS2j  and  S3,  on  the  normal  force  and  pitching  moment  of  store 
S4,  the  bottom  rear  store.  The  circles  indicate  the  influence  of  the 
rear  two  shoulder  stores,  S5  and  S6,  on  the  normal  force  and  pitching 
moment  of  store  Sx . The  inclusion  of  store  S4  in  the  latter  group 
was  not  feasible  because  Sx  was  sting  mounted. 

The  normal-force  coefficient  increments  shown  in  figure  27 (a)  are 
small  in  accordance  with  predictions  of  theory.  The  theory  is  based  on 
determining  the  interference  flow  field  of  the  stores  by  three-dimensional 
source  distributions.  The  pitching-moment  coefficient  increments  for 
Sx  and  S4  have  opposite  signs.  This  effect  is  explained  by  the  theory. 
Stores  Sl5  S2 , and  S3  produce  a downwash  over  the  length  of  S4  which 
decreases  with  distance  behind  the  nose  of  S4.  Accordingly,  a down  load 
is  induced  on  the  nose  of  S4  which  produces  a negative  moment.  Stores 
S 5 and  S6  produce  a downwash  on  store  Sx  which  is  greatest  at  the  rear 
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and  least  at  the  nose.  Accordingly,  the  down  load  at  the  rear  dominates, 
and  a positive  pitching  moment  results.  If  all  of  the  stores  were  sym- 
metrical about  their  mid-points,  the  theoretical  normal-force  increments 
on  Sx  and  S4  would  be  equal,  and  the  moment  coefficient  increments  would 
be  equal  in  magnitude  and  opposite  in  sign  (considering  S4  to  be  inclu- 
ded in  the  rear  grouping) . 

These  data  also  show  that  no  serious  interference  effect  associated 
with  body  doublets  or  trailing  vortices  occurs  up  to  angles  of  attack  of 
10°.  A change  with  angle  of  attack  of  the  increments  would  be  associated 
with  changes  in  store  doublet  strengths  and  changes  in  trail ing-vortex 
characteristics.  Since  the  effects  of  doublets  fall  off  more  rapidly 
with  distance  than  those  of  sources,  the  source  effects  will  dominate 
the  interference.  More  precise  results  can  be  obtained  by  including  the 
doublet  and  trailing-vortex  effects. 

Experimental  results  analogous  to  those  of  figure  27  are  shown  in 

figure  28  for  the  stores  with  cruciform  empennages.  It  is  of  importance 

to  discuss  the  interference  effects  of  Sx  separately  from  those  of  S4 

since  the  empennage  trailing-vortex  systems  of  Sx  , S2,  and  S3  come  close 

to  S4  while  those  of  S5  and  S6  have  negligible  influence  on  Sx . The 

difference  in  the  force  and  moment  acting  on  Sx  (with  empennage)  with 

and  without  S5  and  S6  present  is  due  to  loadings  induced  separately  on 

its  body  and  its  empennage.  If  the  downwash  induced  along  Sx  by  S5  and 

Se  is  not  influenced  significantly  by  their  empennages,  then  the  Sx 

body  loading  contributions  due  to  S5  and  Se  will  be  the  same  with  or 

without  empennages  on  all  stores.  The  rest  of  the  contribution  will  be 

due  solely  to  the  addition  of  the  Sx  empennage  into  the  induced  down- 

wash  field.  Accordingly,  if  the  contribution  in  figure  27  for  Sx  is 

subtracted  from  the  corresponding  body  contribution  in  figure  28,  the 

result  should  be  the  empennage  contribution  of  Sx  due  to  the  addition 

of  S5  and  S6.  When  increments  AC„  and  AC  are  calculated  in  this 

N m 

fashion,  it  is  found  that  the  incremental  center  of  pressure  given  by 
AC^/AC^  lie  on  the  empennage  of  Sx , thus  supporting  the  above  explana- 
tion of  the  interference  effects  of  S5  and  S6  on  S1. 

In  the  assessment  of  the  effects  of  Sx , S2 , and  S3  on  S4,  it  is 
noted  that  the  wakes  of  Sly  S2,  and  S3  all  pass  downstream  along  the 
entire  length  of  S4  so  that  the  downwash  distribution  along  S4  is 
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definitely  influenced  by  the  empennages  of  Sx  , S2  , and  S3.  Accordingly, 
the  center  of  pressure  of  the  interference  effects  determined  in  the 
preceding  manner  will  not  necessarily  lie  on  the  empennage.  In  fact,  it 
is  found  to  lie  approximately  1.5  diameters  behind  the  center  of  moments 
of  S4,  the  store  mid-point. 

7.4  Empennage  Contributions  to  Forces  and  Moments 

Comparison  between  experiment  and  theory  for  forces  and  moments  due 
to  the  presence  of  the  empennage  are  now  considered.  It  is  clear  that 
the  method  of  determining  flow  fields  acting  on  the  store  is  equally 
applicable  to  the  store  empennage  forces  in  a nonuniform  flow  field. 

A method  of  determining  empennage  contributions  to  C^,  CY,  C^,  C^,  and 
C^  based  on  reversed  flow  theorems  has  been  described  in  section  5.3. 
This  method  has  been  used  with  experimental  flow-field  data  to  estimate 
the  above  quantities  for  comparison  with  measurements. 

In  figure  29,  the  forces  and  moments  for  store  Sc  mounted  at  or 
near  the  attached  position,  shown  in  figure  13,  are  shown  as  a function  of 
roll  angle  and  compared  with  theory  for  angles  of  attack  of  8°  and  16°. 

A roll  angle  of  0°  places  the  fins  vertical  and  horizontal.  The  incre- 
mental coefficients  are  the  difference  between  the  coefficients  with  the 
empennage  on  and  with  the  empennage  off. 

The  variations  of  normal-force  and  side-force  coefficient  incre- 
ments with  roll  angle  are  shown  in  figure  29(a) . The  theory  is  in  fair 
agreement  with  the  data  at  both  angles  of  attack.  A measure  of  the 
accuracy  of  the  theory  can  be  made  by  comparing  the  data  at  0=0° 
and  0 = 90°.  In  the  side-force  coefficient  increments,  the  difference 
between  theory  and  experiment  generally  does  not  exceed  the  difference 
between  the  0=0°  and  0 = 90°  experimental  values  of  these  increments. 
Any  difference  can  be  due  to  measurement  inaccuracy  or  to  slight  fin 
misalignment.  We  can  conclude  that  the  agreement  is  within  the  accuracy 
of  the  data. 

With  regard  to  the  comparisons  for  pitching  moment  shown  in 
figure  29(b)  the  data  agrees  with  the  theory  quite  well  when  accounts  are 
taken  of  the  suppressed  origin  of  the  vertical  scale  and  the  difference 
in  the  data  at  0=0°  and  0 = 90°.  With  regard  to  yawing  moment,  the 
difference  between  theory  and  experiment  at  a = 8°  cannot  be  entirely 
accounted  for  by  data  inaccuracy.  At  a = 16°  the  agreement  between 
experiment  and  theory  is  within  the  accuracy  of  the  data. 
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Comparison  between  theory  and  experiment  are  shown  for  rolling- 
moment  coefficient  in  figure  29(c) . In  this  case,  all  the  rolling  moment 
is  due  to  the  empennage.  For  a = 8°  the  comparison  between  experiment 
and  theory  is  fairly  good.  The  difference,  which  at  some  points  amounts 
to  an  increment  in  rolling-moment  coefficient  of  0.02,  corresponds  to 
about  0.35  degrees  of  roll  control  for  a single  opposing  pair  of  panels. 
The  agreement  at  a * 16°  is  not  as  good  with  the  discrepancies  being 
about  twice  those  for  the  a = 8°  case.  The  nature  of  rolling  moment 
for  a cruciform  empennage  is  of  interest  in  assessing  the  comparisons 
between  experiment  and  theory.  In  general,  the  rolling  moment  due  to 
one  opposing  pair  of  panels  is  large  and  opposes  that  due  to  the  other 
pair  of  opposing  panels  with  a small  net  rolling  moment.  Since  the  net 
rolling-moment  coefficients  shown  in  figure  29(c)  are  small,  the  degree 
of  agreement  between  experiment  and  theory  is  considered  fair. 

7.5  Effect  of  Store  Boattail  on  Loadings  and  Forces 

In  the  comparisons  shown  so  far  for  the  single  store,  the  store 
afterbody  has  been  cylindrical  and  this  store  is  designated  Sc*  In 
order  to  study  the  effects  of  boattail  on  the  store  loadings  and  forces, 
a special  pressure  distribution  store  with  boattail,  S^,  as  described  in 
figure  12(d),  but  without  empennage  was  constructed  and  fitted  with 
sufficient  pressure  orifices  to  obtain  accurate  loading  distributions. 

The  nose  of  this  store  is  identical  to  the  nose  of  Sc.  Before  studying 
the  effects  of  boattailing  in  a nonuniform  flow,  it  is  worthwhile  studying 
the  effects  in  uniform  flow. 

7.5.1  Effect  of  store  boattail  on  store  loading  in  uniform  flow 

In  figure  30  the  normal-force  loading  distributions  for  uniform  flow 
for  stores  and  Sc  are  compared  for  several  angles  of  attack  to  show 

the  experimental  effect  of  the  store  boattail  on  the  loading.  Figure 
30(a)  through  30(d)  both  show  identical  experimental  loadings  on  the 
identical  noses,  a fact  indicating  little  or  no  upstream  influence  of 
the  separated  flow  regions  downstream  of  the  nose.  However,  the  boat- 
tailed  afterbody  carries  more  negative  load  than  the  cylindrical  after- 
body as  might  be  expected. 

Several  theories  are  shown  on  the  figure  for  comparison  with  the 
data.  Slender-body  theory,  section  5.2.2,  does  a fair  job  of  calculating 
the  load  distribution  at  a = 4°.  In  order  to  obtain  a more  precise 
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theoretical  prediction  for  attached  flow  for  assessing  separation  effects, 
a fully  three-dimensional  source  and  doublet  calculation,  D.  S.,  was 
made  by  the  method  of  Appendix  III.  The  question  of  the  Kutta  condition 
for  a body  of  revolution  was  ignored,  and  the  blunt  base  of  the  body  was 
continued  smoothly  to  a point.  It  is  presumed  that  this  solution  will 
give  an  accurate  estimate  of  the  load  distribution  for  the  body  except 
possibly  close  to  the  base.  The  boundary  conditions  of  thickness  and 
angle  of  attack  were  satisfied  at  the  body  surface,  and  the  full  Bernoulli 
equation  was  used  in  determining  the  loadings. 

It  can  be  seen  that  the  data  are  in  excellent  agreement  with  the 
doublet  theory  up  to  the  body  separation  region.  Assuming  that  the 
separation  region  for  the  present  purposes  can  be  taken  as  that  region 
where  the  experimental  loading  departs  from  the  theoretical  one,  one 
finds  a smooth  forward  movement  of  the  separation  region  with  increase 
in  angle  of  attack  as  shown  in  figure  31(a).  In  this  figure  the  separa- 
tion regions  are  shown  as  estimated  from  figure  30.  A mean  line  has  been 
faired  through  the  regions.  The  mean  value  of  xg  from  this  curve  has 
been  used  in  determining  the  curve  in  figure  31(b),  which  shows  the 
distance  of  the  separation  point  behind  the  shoulder  of  the  store  in 
multiples  of  the  maximum  diameter.  Separation  lies  within  about  1 to  2 
maximum  diameters  behind  the  store  shoulder  for  the  boattail  body  for 
most  of  the  useful  angle-of-attack  range 

The  source-doublet  theory  is  also  given  in  figure  30  for  the  store 
with  the  cylindrical  afterbody.  The  separation  parameters  have  been 
determined  from  these  data  and  are  shown  on  figure  31  for  comparison  with 
those  for  the  boattail  store.  It  is  noted  that  boattail  moves  the 
separation  region  forward  as  would  be  expected. 

With  regard  to  the  additional  loading  behind  the  separation  point 
due  to  cross-flow  drag,  analysis  of  the  data  does  not  correspond  to  a 
uniform  cross-flow  drag  coefficient.  It  has  been  found  that  an  estimate 
of  body  loading  can  be  made  by  assuming  the  loading  to  be  uniform 
downstream  of  separation,  at  its  separation  value.  A slight  effect  of 
the  base  is  evident  very  close  to  the  base. 
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7.5.2  Effect  of  boattail  on  store  loading  in  presence  of  wing-body 
combination 

Normal-force  and  side-force  loadings  were  measured  on  the  store  in 
the  carriage  position  in  the  presence  of  the  wing-body  combination,  see 
figure  13,  without  the  pylon  present.  The  normal-force  loading  distri- 
butions with  and  without  boattail  are  shown  in  figure  32  for  a = 4° 
and  a = 16°.  The  theory  shown  for  a = 4°  is  the  slender-body  theory 
which  predicts  quite  well  the  loading  distribution  with  and  without 
boattail.  The  theory  is  based  on  theoretical  flow  angles  at  a = 4° 
so  that  the  agreement  between  experiment  and  theory  includes  the  errors 
in  the  flow-field  prediction  method  and  the  loading  calculation  method. 

At  a « 16°  the  data  as  shown  in  figure  32(b)  exhibit  a small 
effect  of  boattail.  On  the  nose  the  results  for  both  stores  would  be 
identical  except  for  possible  upstream  influence  of  the  afterbody  or 
possible  experimental  error.  The  experiments  were  made  many  months  apart. 
It  is  noted  that  slender-body  theory  over-predicts  the  magnitude  of  the 
loading  distribution  at  this  large  angle  of  attack. 

The  differences  between  the  data  and  slender-body  theory  are  due 
in  part  to  an  overprediction  of  the  angle-of-attack  distribution  along  the 
store  and  in  part  to  an  overprediction  in  the  loading  distribution  calcu- 
lation. To  eliminate  both  these  sources  of  error  in  the  theory  a source- 
doublet  calculation,  Appendix  III,  was  made  similar  to  that  for  uniform 
flow  previously  described  in  connection  with  figure  30  using  the  experi- 
mental downwash  distribution  as  input.  The  improved  theory  indicates 
generally  improved  agreement  with  experiment  especially  when  the  buoyancy 
correction  is  included  in  the  theory.  This  was  made  using  the  method  of 
section  5.2.1.  These  comparisons  are  of  interest  because  they  tend  to 
show  that  separation  seen  for  uniform  flow  is  of  reduced  extent  for  the 
store  at  the  same  angle  of  attack  under  the  wing-body  combination.  This 
result  can  be  explained  by  the  fact  that  the  wing-body  combination  turns 
the  flow  more  parallel  to  the  store  and  reduces  the  average  angle  of 
attack  in  the  attached  position.  For  instance,  with  a = 16°,  the  local 
angle  of  attack  at  the  nose  is  about  13°  while  that  on  the  rear  half  of 
the  store  is  about  6°.  Whether  there  is  any  separation  over  the  rear  part 
of  the  store  cannot  be  determined  with  certainty  because  the  accuracy  of 
the  improved  theory,  Appendix  III,  for  attached  flow  cannot  be  guaranteed 
in  the  region  of  the  body  base  because  of  the  unknown  wake  effects  and  the 
lack  of  a Kutta  condition  applicable  to  bodies  of  revolution. 
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Once  the  store  leaves  the  attached  position,  it  is  not  clear  that 
the  angle-of-attack  distribution  along  the  store  will  be  small.  Large 
angles  of  attack  may  occur  once  the  straightening  effect  of  the  wing  is 
reduced,  and  extensive  separation  similar  to  that  in  uniform  flow  may 
occur.  In  this  case,  a method  of  calculating  the  loading  distribution 
similar  to  that  suggested  for  uniform  flow  may  be  worthwhile.  This  point 
has  not  been  pursued  in  the  current  investigation  since  an  accurate  method 
of  determining  the  wing  flow  field  at  high  angles  of  attack  is  unavailable. 
Data  are  available  for  a boattail  store  at  distances  up  to  a diameter 
below  the  carriage  position  but  no  further. 

Results  similar  to  those  in  figure  32  for  normal  force  are  shown 
in  figure  33  for  side  force.  At  a = 4°  the  effects  of  boattail  on 
side-force  loading  distribution  are  small  and  are  well  predicted  by 
s lender-body  theory.  At  a = 16°  slender-body  theory  overpredicts  the 
magnitudes  of  the  loading  distribution.  The  use  of  the  more  precise  doublet- 
source  method  yields  better  agreement  with  experiment  although  the  correction 
for  buoyancy  tends  not  to  improve  the  agreement.  The  sidewash  angle 
varies  from  9.2°  at  the  nose  to  4.3°  at  the  base  in  a nearly  linear 
fashion.  This  decreasing  local  sidewash  angle  in  the  strearnwise  direction 
probably  tends  to  inhibit  separation  although  the  extent  of  separation, 
if  any,  cannot  be  determined  precisely  from  the  results  of  the  figure. 

7.5.3  Effect  of  boattail  on  empennage  forces  and  moments 

All  forces  and  moments  except  drag  have  been  measured  on  and 

with  and  without  empennages.  The  changes  in  C^,  CY,  C^,  and  due 

to  adding  the  empennage  are  shown  in  figure  34  as  a function  of  roll 
orientation.  In  these  tests  the  pylon  was  present,  and  the  store  is 
slightly  below  the  pylon,  Az/d  = 0.1333.  The  experimental  differences 
due  to  boattail  can  be  seen  by  comparing  solid  and  open  symbols.  The 
theory  is  based  on  measured  flow  quantities  in  the  region  to  be  occupied 
by  the  empennage  together  with  reverse  flow  theorems. 

For  normal  force  and  side  force,  the  effect  of  boattail  is  to 
increase  the  magnitude  of  the  coefficients  by  an  amount  in  good  accordance 
with  theory.  The  discrepancy  between  experiment  and  theory  is  within 
the  accuracy  of  the  experimental  data  judged  by  comparing  values  for 
0 = 0°  and  <J>  = 90°,  which  should  be  equal.  The  increase  in  the  magnitude 
of  the  normal  force  due  to  boattail  is  about  12  percent  and  about  30  per- 
cent for  side  force. 
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The  results  for  pitching  moment  and  yawing  moment  shown  in  figures 
34(c)  and  34(d)  are  in  good  accord  with  theory  when  account  is  taken  of 
the  suppressed  origin  and  the  inaccuracy  of  the  data  exhibited  by  compari- 
son of  the  values  at  <£  = 0°  and  $ = 90°. 

The  results  for  rolling  moment  shown  in  figure  34(e)  exhibit  the 
same  behavior  as  those  in  figure  29(c). 

7.6  Trajectory  Studies 

The  methods  of  flow-field  predictions  and  store  force  and  moment 
determinations  presented  in  sections  4 and  5 have  been  combined  into  one 
program  with  the  six-degree-of-f reedom  equations  of  motion  of  section  6 
to  yield  a trajectory  prediction  method.  A number  of  sample  trajectories 
have  been  run  to  provide  predictions  for  comparison  with  captive-store 
trajectories . 

In  order  to  represent  full-scale  conditions,  the  wind-tunnel  models 
shown  in  figures  7,  11,  and  12  have  been  scaled  up  by  a factor  of  twenty. 
The  following  input  quantities  have  been  used  for  all  trajectories  to  be 
shown: 

Store  Mass  = 15.53  slugs 
I =8  slug-ft2 

XX  ^ 


I = I =80  slug-ft2 
yy  zz 


I = I = I =0 
xy  yz  xz 


Altitude  = 5000  ft 

£ = 10  ft/sec  (initial  vertical 

velocity  of  ejected  store) 

C t=  = 0.625  ft  (initially,  ejected 
store  is  one  radius  below 
attached  position) 

Initial  store  and  aircraft  angles  of  attack  are  equal. 
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7.6.1  Statically  unstable  single  store  dropped  beneath  wing 

A comparison  between  the  calculated  trajectory  and  the  captive- 
store  trajectory  for  no  damping  is  shown  in  figure  35  for  the  cylindrical 
afterbody  store  (fig.  12(c))  with  no  empennage.  It  is  released  at  time 
t = 0 one  store  radius  beneath  its  attached  position  on  the  pylon  at  the 
one-third  semispan  position,  figure  13,  with  a = 4°  for  both  the  air- 
craft and  the  store.  In  figure  35  the  three  left-hand  curves  show  the 
position  of  the  store  center  of  gravity  relative  to  its  carriage  position. 
The  quantities  A£ , Ar),  and  A £ are  respectively  positive  forward, 
positive  to  the  right  (toward  the  fuselage) , and  positive  down.  The  left- 
hand  curves  show  a slight  rearward  movement  of  the  store,  no  lateral 
movement,  and  a vertical  movement  equivalent  to  free  fall.  The  two  right- 
hand  curves  show  a substantial  yaw  AY  of  the  store  nose  outboard 
and  a large  pitch  up  A0  of  the  store  which  is  somewhat  overpredicted 
by  the  computer  program. 

The  fact  that  the  pitch  up  is  overpredicted  can  be  readily  explained. 
The  nose-up  motion  of  the  store  places  it  at  an  angle  of  attack  of  about 
12°  at  t = 0.4  seconds.  At  this  time  the  store  is  well  below  the  wing 
so  that  the  flow  is  nearly  uniform.  At  this  condition  the  captive  store 
balance  yields  a pitching-moment  coefficient  of  0.9,  the  store  alone  in  a 
uniform  field  yields  1.05,  and  the  computer  program  yields  1.3.  It 
appears  that  the  captive-store  measurements  are  too  low  and  the  computer 
program  result  is  too  high.  Figure  30(c)  shows  the  store  loading  at 
a = 12°  in  uniform  flow.  It  is  seen  that  slender-body  theory  overpredicts 
the  loading  on  the  store  nose  at  this  angle  of  attack  and  that  the  contri- 
bution of  the  afterbody  to  the  pitching  moment  is  small.  Increased  accuracy 
of  prediction  would  then  be  obtained  using  the  improved  source-doublet 
method  described  in  connection  with  figure  30. 

The  effect  of  damping  for  this  statically  unstable  case  is  of 
little  interest  and  was  not  studied.  The  calculated  trajectory  shown 
required  approximately  5 to  6 minutes  on  a Univac  1108  computer. 

7.6.2  Statically  stable  single  store  dropped  beneath  wing 

A comparison  similar  to  that  of  figure  35  is  shown  in  figure  36  for 
the  same  store  and  initial  conditions  except  that  a cruciform  empennage 
has  been  added  to  produce  static  stability.  Results  are  shown  for  angles 
of  attack  of  a = 0°  and  a = 4°.  The  store  center  of  gravity  coordinates 
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(At;,  Ar|,  AC)  show  the  same  general  behavior  as  the  previous  case.  The 
pitching  and  yawing  motions  from  the  computer  program  and  the  captive- 
store  trajectory  are  in  good  agreement  for  both  angles  of  attack.  The 
pitching  motion  is  nose  down  in  this  instance  compared  to  a pitch-up  in 
the  previous  case.  It  is  noted  that  a few  degrees  of  roll  occur.  The 
roll  excursions  are  so  small  that  they  have  no  effect  on  the  calculated 
trajectories . 

In  order  to  assess  the  importance  of  damping  on  the  motion,  the 
trajectory  was  also  calculated  including  damping  in  all  three  angular 
motions.  The  damping  effect  is  negligible  except  in  the  pitching 
oscillation  which  was  reduced  in  maximum  amplitude  by  about  one  degree 
as  the  result  of  pitch  damping. 

7.6.3  Statically  unstable  store  dropped  from  TER  group  beneath 
wing 

Stores  Sc  without  empennages  were  mounted  in  a TER  group  under 
the  wing,  figure  18,  and  the  trajectory  of  the  bottom  store  was  calcu- 
lated using  the  computer  program  and  measured  by  the  captive-store 
technique.  The  comparison  between  the  two  sets  of  results  is  shown  in 
figure  37.  As  in  the  preceding  case,  the  center-of-gravity  location 
shows  the  same  behavior  with  the  same  good  agreement  between  theory  and 
measurement.  The  store  exhibits  divergence  in  both  pitch  and  yaw.  The 
overprediction  of  the  pitch  up  in  this  case  is  similar  to  that  shown  in 
figure  35  for  the  same  reasons. 

8.  CONCLUDING  REMARKS 

This  report  presents  the  results  of  an  investigation  which  has  been 
conducted  with  the  objectives  of  extending  and  improving  the  store 
separation  prediction  method  of  reference  1.  In  this  earlier  work  the 
trajectory  was  limited  to  three  degrees  of  freedom,  the  pylon  and  ejector 
rack  were  not  modeled,  wing  dihedral  was  not  included,  a two-dimensional 
wing  thickness  model  was  utilized,  and  a simple  method  was  used  to  calcu- 
late the  empennage  normal  force  which  could  not  be  used  to  determine 
rolling  moment.  All  of  these  shortcomings  have  been  eliminated  in  the 
present  work. 

The  fuselage  and  store  volumes  are  represented  by  distributions  of 
three-dimensional  point  sources  along  the  body  longitudinal  axes.  The 
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same  method  is  used  to  model  the  body  of  the  ejector  rack.  The  short 
pylons  on  the  rack  have  not  been  included  in  the  present  method. 

The  wing  model  used  in  the  work  of  reference  1 has  been  replaced  by 
a wing-pylon  model  which  includes  complete  mutual  interference  between 
the  wing  and  pylon.  The  wing-pylon  thickness  distribution  is  represented 
by  three-dimensional  source  strips  and  superimposed  upon  this  is  a 
vortex  lattice  which  represents  the  loading  distribution.  The  boundary 
condition  imposed  in  determining  the  vortex  strengths  includes  wing 
angle  of  attack  and  twist  and  camber,  mutual  interference  between  wing 
and  pylon,  and  interference  due  to  the  fuselage,  rack,  and  stores.  The 
wing  model  allows  for  breaks  in  sweep  and  dihedral.  The  pylon  can  be 
swept  but  breaks  in  sweep  are  not  allowed. 

Calculations  utilizing  the  flow  models  just  described  have  been 
compared  with  experimental  flow-field  data  obtained  under  various  parent 
aircraft  configurations.  The  agreement  is  quite  good  except  for  the  case 
of  a store  in  the  carriage  position  on  a TER  rack  where  the  short  pylon 
which  is  part  of  the  rack  produces  an  effect  not  accounted  for  by  the 
theory.  This  effect  disappears  rapidly  as  the  store  moves  off  of  the 
rack.  This  pylon  can  be  accounted  for  within  the  framework  of  the  present 
theory  but  has  not  at  this  time. 

Methods  are  presented  for  calculating  all  of  the  aerodynamic  forces 
and  moments  except  axial  force.  The  methods  for  calculating  the  body 
forces  and  moments  are  the  same  as  those  used  in  reference  1 except  that 
a simpler  method  based  on  slender-body  theory  is  used  for  the  buoyancy  forces 
and  moments.  Angle  of  attack  effects  are  accounted  for  by  a combination 
of  slender-body  theory  and  viscous  crossflow  theory.  The  empennage  forces 
and  moments,  including  rolling  moment,  are  calculated  using  reverse  flow 
theorem.  Good  agreement  between  calculated  and  experimental  forces  and 
moments  is  obtained  except  at  high  angles  of  attack  where  slender  body 
theory  overpredicts  the  lift  on  the  store  nose.  It  is  shown  that  by 
using  a more  precise  method  for  calculating  the  forces  the  agreement  is 
significantly  improved.  For  a group  of  six  stores  mounted  on  a MER  rack 
the  theory  predicts  the  interference  of  one  cluster  of  three  on  the  other 
cluster  quite  well. 

Theory  and  experiment  have  been  used  to  determine  the  effect  of 
store  boattail  on  the  location  of  flow  separation  from  the  body. 
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The  results  indicate  that  the  effect  of  boattail  is  to  cause  a forward 
movement  of  the  separation  location  for  the  store  in  uniform  flow.  With 
a parent  aircraft  the  effect  is  not  as  pronounced  since  the  wing  tends  to 
align  the  flow  with  the  body  axis  thus  reducing  the  average  angle  of 
attack  over  the  rear  half  of  the  store. 

The  flow  field  and  force  and  moment  prediction  methods  have  been 
included  with  the  six-degree-of-freedom  equations  of  motion  in  a computer 
program  which  calculates  the  trajectory  of  a store  dropped  from  an  air- 
craft flying  at  constant  velocity  in  rectilinear  flight.  The  stores  can 
be  mounted  singly  on  a pylon  or  in  groups  on  a rack.  The  equations  of 
motion  allow  for  mass  and  inertia  assymmetries  of  the  ejected  store. 
Comparisons  have  been  made  between  predicted  trajectories  and  captive- 
store  trajectories  obtained  in  the  wind  tunnel.  Very  good  agreement  is 
obtained  for  stable  stores.  For  unstable  stores  very  good  agreement  is 
shown  except  for  the  store  pitch  angle.  The  theory  overpredicts  the 
pitch  excursion  from  its  initial  value  due  to  slender-body  theory  over- 
predicting the  load  on  the  store  nose. 
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Figure  1.-  Coordinate  system  for  an 
axisymmetric  body. 
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Figure  2.-  Coordinate  systems  for  wing-pylon 
and  horseshoe  vortices. 
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Figure  3.-  Vortex-lattice  model  of  a wing-pylon. 


Figure  4 - Flow  tangency  boundary  conditions 
for  the  wing  and  pylon. 
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Figure  5.—  Coordinate  system  for  wing  thickness  source  strips. 


62 


Figure  6.-  Coordinate  system  for  pylon  thickness  source  strips. 
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Figure  7.-  Details  of  TER  rack. 
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Figure  8.-  Coordinate  «ystems  fixed  in  ejected  store  and  used 
in  force  and  moment  calculation. 
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Figure  9.-  Coordinate  systems  used  in  trajectory  calculation. 
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Figure  10.-  Coordinate  systems  used  in  empennage  force 

and  moment  calculation. 
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Figure  11.-  Wing-fuselage  combination. 
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For  wing  pylons,  pylon  centerline  located  at 
40%  wing  chord. 

For  fuselage  pylon,  pylon  centerline  located 
19.43  inches  aft  of  fuselage  nose. 


(a)  Pylon  used  with  single  store  and 
MER  and  TER  racks. 


Figure  12.-  Aircraft  components  used  with 
wing- fuselage  combination. 
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(b)  MER  rack. 


Figure  12.-  Continued. 


(c)  Large  store  with  cylindrical  afterbody 
and  rectangular  cruciform  fins. 

Figure  12.-  Continued. 


(d)  Large  store  with  boattail  afterbody 
and  rectangular  cruciform  fins. 


Figure  12.-  Continued. 
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All  dimensions 
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Store  midpoint  aligned  with  suspension  centerline 
on  MER  rack  when  in  carriage  position. 

Unfinned  force  and  moment  model  same  as  finned 
model  except  for  omission  of  fins. 

Finned  and  unfinned  dummy  stores  same  as  respective 
force  and  moment  stores  except  for  omission  of 
center  bore  for  balance  mounting. 


(e)  Small  boattail  MER  store  with 
cruciform  empennage. 


Figure  12.-  Concluded. 
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Figure  13.-  Large  store  with  cylindrical  afterbody 
under  wing-fuselage  configuration. 
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Figure  14.-  Effect  of  pylon  on  flow  field  of 
wing-body  combination  at  centerline 
location  of  attached  store. 
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Figure  15. - Effect  of  pylon  on  load  distributions 
of  attached  store. 
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Pylon  • 


Figure  15.-  Concluded. 
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Figure  16.-  Effect  of  vertical  location  on  normal-force 
and  side-force  coefficients  of  store. 
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Figure  18.-  Large  stores  with  cylindrical  afterbody 
in  a TER  arrangement  under  the  wing. 
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Figure  19.-  Effect  of  vertical  position  on  downwash 
and  sidewash  under  TER  configuration 
(parent  aircraft  - WBPTS2S3) . 
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Figure  20.-  Effect  of  angle  of  attack  on  downwash 
and  sidewash  under  TER  configuration 
(parent  aircraft  - WBPTS2S3) . 
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(a)  Normal-force  distribution. 

Figure  21.-  Effect  of  vertical  position  on  loading  distributions 
on  store  Sj  of  TER  configuration 
(Parent  aircraft  WBPTS2S3) . 
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(b)  Side-force  distribution. 
Figure  21.-  Concluded. 
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a = 6°  • 


(a)  Normal-force  distribution. 


Figure  22.-  Effect  of  angle  of  attack  on  loading  distributions 
on  Si  in  TER  configuration  (parent  aircraft  - WBPTS2S3) . 
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(b)  Side-force  distribution. 
Figure  22.-  Concluded. 
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Without  TER  • 


(a)  Normal-force  distribution* 

Figure  23.-  Effect  of  TER  on  loading  of  store  Sx 
(parent  aircraft  - WBPT  and  WBP) . 
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Figure  23.-  Concluded. 
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Figure  24.-  Effect  of  vertical  position  on  loading 
distribution  of  shoulder  store  S2 
of  TER  configuration 
(Parent  aircraft  WBPTS3) . 
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Figure  24.-  Concluded. 
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a = 6°  • 


(a)  Normal-force  distribution. 

Figure  25.-  Effect  of  angle  of  attack  on  loading  distribution 
of  shoulder  store  S2  of  TER  configuration 
(parent  aircraft  - WBPTS3) . 
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a = 6°  • 


Figure  25.-  Concluded. 
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Figure  26.-  Small  boattail  stores  in  a MER  arrangement 

under  the  fuselage. 
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STORE  DATA  THEORY 


Figure  27.-  Effects  of  S5,  Se  on  Sx  and  Sx , S2,  S3 
on  S4  in  MER  configuration;  no  empennages. 
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(b)  Pitching-moment  coefficient 

Figure  28.-  Effects  of  S5,  S6  on  Sx  and  Sl9  S29  S3 
on  S4  in  MER  configuration;  with  empennages. 
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(a)  Normal-force  and  side-force  coefficients. 


Figure  29.-  Effect  of  roll  angle  on  empennage  contributions 

to  store  S forces  and  moments 
c 

in  the  attached  position 
(parent  aircraft  WBP) . 
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a = 16°;  Az/d  = 0.1333 
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(b)  Pitching-moment  and  yawing-moment  coefficients. 
Figure  29.-  Continued. 
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Figure  29.-  Concluded. 
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S.B.  - SLENDER  BODY 
D.S.  - DOUBLET  SOURCE 


Figure  30.-  Effects  of  store  boattail  on 
store  loading  in  uniform  flow. 
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S.B.  - SLENDER  BODY 
D.S.  - DOUBLET  SOURCE 


(b)  a = 8°. 

Figure  30.-  Continued. 
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S.B. -SLENDER  BODY 


(c)  a = 12°. 
Figure  30.-  Continued. 
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(d)  a - 16°. 
Figure  30.-  Concluded. 
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(a)  Separation  location. 


Figure  31.-  Effects  of  boattail  on  location 
of  separation  on  store  afterbody; 
uniform  flow. 
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Figure  32.-  Effects  of  boattail  on  normal-force 
distribution  of  store  in  presence 
of  wing-body  combination. 
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Figure  32.-  Concluded. 
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(a)  a = 4°. 

Figure  33.-  Effects  of  boattail  on  side-force 
distribution  of  store  in  presence 
of  wing-body  combination. 
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Figure  33.-  Concluded. 
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(a)  Normal-force  coefficient. 

Figure  34.-  Effects  of  boattail  on  the  forces 
and  moments  due  to  the  empennage 
(parent  aircraft  WBP) . 
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Figure  34.-  Continued. 
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(c)  Pitching-moment  coefficient. 
Figure  34.-  Continued. 


110 


DATA  THEORY 


Sc  o 


Figure  34.-  Continued. 
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(e)  Rolling-moment  coefficient. 
Figure  34.-  Concluded. 
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CAPTIVE  STORE  THEORY 
NO  DAMPING  NO  DAMPING 


af  * 4°  o 


Figure  35.-  Comparison  between  calculated  trajectory  and  captive-store 
trajectory  of  store  S«  released  at  one-third 
semispan  location;  - 0.4. 
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NO  DAMPING  THEORY 

af  * O*  O NO  DAMPING 

af  - 4*  □ NO  DAMPING 

af  « 4°  WITH  DAMPING 


Figure  36.-  Comparison  between  calculated  trajectories  and  captive-store 
trajectories  of  store  Sc  with  empennage  released  at 
one-third  semispan  location;  Hm  - 0.4. 
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CAPTIVE  STORE  THEORY 
NO  DAMPING  NO  DAMPING 


af  « 4°  o 


Figure  37.-  Comparison  between  calculated  trajectory  and  captive-store  trajectory 
of  finless  store  Sj  released  from  TER  group  at  one-third 
semispan  location;  Mm  - 0.4. 
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APPENDIX  I 


CALCULATION  OF  FORCES  AND  MOMENTS  DUE  TO  THE 
EMPENNAGE  OF  AN  EXTERNAL  STORE  IN  A NONUNIFORM  AIRFLOW 

1-1.  INTRODUCTION 

The  empennage  of  an  external  store  in  a nonuniform  airstream,  as 
may  occur  under  an  aircraft,  will  develop  forces  and  moments  which  affect 
its  trajectory.  While  the  forces  on  the  body  of  the  store  can  be  deter- 
mined using  slender-body  theory,  the  empennage  is  generally  not  slender  so 
that  a different  approach  is  required.  The  approach  here  will  be  based  on 
a combination  of  reverse  flow  theorems  and  slender  body  theory  scaled  for 
the  effects  of  aspect  ratio  as  done  in  wing-body  interference  theory 
(ref.  16) . It  will  be  possible  to  calculate  all  forces  and  moments 
except  drag  using  a unified  approach.  The  effects  of  interference  between 
body  and  tail  panels  will  be  accounted  for,  and  the  method  will  be 
applicable  to  cases  where  vortices  appear  in  the  flow;  for  instance, 
vortices  discharged  from  the  body  at  a position  in  front  of  the  empennage. 
Also  planar,  triform,  and  cruciform  empennages  are  included  within  the 
framework  of  the  method.  However,  the  influence  functions  for  triform 
empennages  will  not  be  worked  out  in  this  appendix.  Only  planar  and 
cruciform  empennages  will  be  treated. 

1-2.  AXIS  SYSTEM  AND  NOTATION 

The  coordinate  system  to  be  used  in  the  empennage  calculation  is 
shown  in  figure  10  for  the  crossflow  plane  at  which  the  empennage  forces 
are  assumed  to  act,  the  quarter-chord  position  of  the  mean  aerodynamic 
chord  of  the  exposed  panels.  This  xf»Vf>zf  system  is  fixed  in  the  fins 
and  oriented  with  respect  to  the  store  body  coordinate  system  of  figure  8 
through  the  angle  0f.  The  xf  axis  lies  on  the  xg  axis.  The  angle 
is  specified  in  order  that  the  tail  fins  can  be  other  than  vertical 
and  horizontal  at  the  beginning  of  a trajectory.  Both  coordinate  systems 
are  fixed  in  the  store  so  that  they  translate  and  rotate  with  the  store 
during  the  trajectory. 

The  cruciform  fins  are  designated  Fly  F2 , F3,  and  F4  as  shown  in 
figure  10  with  F3  and  F4  missing  for  the  planar  case.  The  resultant 
velocity  components  normal  to  the  four  fins  are  designated  Wly  W2 , V1 , 
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and  V2  , respectively,  with  positive  directions  as  shown.  They  vary 
across  the  span  of  each  fin.  The  procedure  for  calculating  these  normal 
velocities  is  as  follows: 

1.  Specify  panel  coordinates  (y^  or  z^)  where  the  normal  velocities 
are  to  be  obtained. 

2.  Transform  there  points  to  the  store  body  coordinate  system  of 
figure  8. 

3.  Calculate  the  Vg  and  Wg  velocities  by  the  procedure  described 
in  section  5.1. 

4.  From  Vg  and  for  each  point  determine  the  component  normal 

to  the  fin. 

The  velocity  determined  in  the  above  manner  include  terms  accounting 
for  the  store  translational  motion  and  pitch  and  yaw  rotation.  The  method 
of  accounting  for  rolling  motion  will  be  discussed  in  section  1-4  of  this 
appendix . 


1-3.  DETERMINATION  OF  NORMAL  FORCE  AND  SIDE  FORCE  DUE  TO  EMPENNAGE 

The  normal  force,  Z^,  and  side  force,  Y^,  due  to  the  nonuniform 
flow  field  at  the  empennage  is  readily  determined  from  reverse  flow 
theorems.  The  calculation  is  identical  whether  the  empennage  is  a 
planar  one  with  panels  lying  on  the  y^  axis,  or  a cruciform  one  with 
panels  on  both  the  y^  and  z^  axes.  In  determining  Z^  for  the 
planar  case,  we  assume  that  the  velocities  Wx  and  W2  depend  principally 
on  y^  and  that  they  vary  slowly  enough  along  the  chordwise  direction 
that  they  may  be  represented  by  an  average  value.  Wind-tunnel  measure- 
ments bear  out  this  assumption.  The  values  of  Wx  and  W2  yield  symmetri- 
cal and  unsymmetrical  angles  of  attack  as  follows: 


as  = 


wi  <yf)  + w2  (-yf) 


2V 


a = 
u 


wi (yf)  - w2 (-yf) 


2V 


j 


a-i) 


where  is  the  velocity  of  the  store  moment  center  which  is  assumed 

s 
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equal  to  the  resultant  velocity  at  the  points  in  question.  The  symmetri- 
cal part  only  produces  Z^,  and  the  unsymmetrical  part  produces  only 
rolling  moment.  By  symmetry  there  is  no  effect  of  panels  F3  and  F4 
in  the  production  of  normal  force  by  panels  F1  and  F2 . However  panel- 
panel  interference  is  significant  for  rolling  moment.  The  lack  of  panel- 
panel  interference  for  and  Zf  makes  the  calculation  of  these 

forces  identical  for  planar  and  cruciform  configuration. 

Consider  now  a planar  empennage  with  a circular  body  in  direct  flow 
and  in  reverse  flow  as  shown. 


Let  S be  the  body  planform  area  and  S be  the  total  planform  area 
of  both  panels.  Let  P^  and  be  the  pressure  coefficient  and  local 

angle  of  attack  in  direct  flow,  and  let  Pr  and  af  be  the  corresponding 
quantities  in  reverse  flow.  Then  by  the  reverse  flow  theorem  (ref.  16, 

p.  221) 


P a ,dS  = 
r d 


P ,a  dS 
d r 

+ SH 


(1-2) 


Choose  the  following  specific  values  of  angles  of  attack  and  pressure 
coefficients 


a , = 0 on  S_  a = 0 on  S_ 
a b r B 


ad  = as  on  SH  ar  = 1 on  SH 


P,  = Pj 
d d 


P._  = P i 
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Then 


(1-3) 


Now  the  normal  force  on  the  horizontal  panels  in  the  presence  of  the 
body  considering  both  upper  and  lower  surfaces  is 


Sh  fc-e- 


,H(B)  = 2c>«s  JJP ddS  = ^ f J as  Pi  dxf  d^f  (I-4) 


H 


a 1 . e . 


Let  the  span-load  distribution  in  reversed  flow  be  (cc J such  that 

* i 


Accordingly  we  have 


(cep  x = 2 


t . e . 

/ p> 


dx. 


l.e. 


sh 

‘ f , H (B)  = / as(CCi)1  d*i 


(1-5) 


(1-6) 


since  as  depends  only  on  y^. 

We  now  establish  the  amount  of  normal  force  carried  over  onto  the 
body  by  viture  of  the  panel  pressure  field.  For  this  case  let 


Then 


a , = 0 on  S a = 1 on  S_ 

d B r B 


ad  = as  °n  SH  ar  = 0 °n  SH 


P . = P, 
d d 


P..  = Pr 


ff  PsasdSH 

SH 


(1-7) 
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The  lift  on  the  body  due  to  the  horizontal  panels  is  now 


■ If* ddSB  ' 4q«s  /h  / “sp2  dXf  dy£ 


(1-8) 


a l.e. 


Designating  the  span  loading  for  this  case  by  (cc J we  have 

^ 2 


so  that 


t.e. 


icc«>2  - 2 


/ 


P2  dxf 


l.e. 


'f.B(H) 


rh 

= 2%,  a 

s J 


dyf 


(1-9) 


(1-10) 


The  total  normal  force  due  to  the  empennage  in  the  z^  direction  is 
designated  (Z0H  - Zfi)  and  is  given  by 


(ZBH  ZB*  f = Zf,H(B)  + Zf,B(H) 


(cc 


l 


f 


- 2Clcc 


/h 


(cci)3  °s 


d-ii) 


where  (cc J is  the  span-loading  on  the  horizontal  panel  for  the  configu- 

Ju  3 

ration  in  reverse  flow  with  both  body  and  panels  at  unit  angle  of  attack. 

In  the  application  of  equation  (1-11)  to  the  empennage,  we  now 
make  an  approximation  which  simplifies  the  calculation  greatly  with 
little  loss  of  accuracy.  If  the  value  of  (cc^)3  from  slender-body  theory 
is  used  in  equation  (1-11) , we  obtain  the  slender-body  approximation  to 
the  normal  force.  Since  (cc J is  the  same  in  direct  and  reverse  flow, 

Ju  3 

it  is  easily  obtained.  However,  we  know  that  slender-body  theory  tends 
to  overpredict  the  normal  force  by  an  amount  which  increases  with  the 
aspect  ratio  of  the  lifting  surface  formed  by  joining  the  exposed  tail 
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panels  together,  the  so-called  tail  alone.  It  has  been  found1  that  we 

can  correct  the  slender  body  result  for  this  aspect  ratio  effect  by 

multiplying  it  by  the  ratio  of  the  lift-curve  slope  of  the  tail  alone  on 

the  basis  of  linear  theory  to  that  on  the  basis  of  slender-body  theory. 

Let  (dC^/daJ H be  the  lift-curve  slope  of  the  horizontal  surfaces 

joined  together  on  the  basis  of  linear  theory  with  an  arbitrary  reference 

area  S_.  Then  the  ratio  of  the  lift-curve  slopes  is  simply 

K ri 


2Tr(sh  -a)2 


(1-12) 


Based  on  the  same  reference  area,  the  normal  force  coefficient  corrected 
for  aspect  ratio  is 


BH 


(Z 


BH 


s 


R 


as(cci) 


dyf 


(1-13) 


where  (cc^)^  is  the  slender-body  theory  span  loading  given  by 


(cci>3 


-4V 


( sh2 *yf2  ~ a4) (sh2  - yf2 ) 


s 2v  2 
n yf 


(1-14) 


from  reference  16,  equation  (6-39) . Equation  (1-13)  is  valid  provided 
Cz  and  ^dCL/do^H  are  kased  on  t'*ie  same  reference  area.  If  ^dCL/da)  H 
is  per  degree,  then  ag  should  be  in  degrees. 

In  application  of  equation  (1-13)  the  velocities  Wx (yf)  and  w2(-yf) 
are  determined  from  flow  field  calculations  without  the  store  present. 
Because  of  the  store  the  values  of  Wx (yf)  and  W2(-yf)  will  be  enhanced  by 
"Beskin  upwash."  Thus,  if  we  assume  the  store  is  subject  to  a crossflow 


1Pitts,  W.  C.,  Nielsen,  J.  N.,  and  Kaattari,  G.  E.:  "Lift  and  Center  of 

Pressure  of  Wing-Body-Tail  Combinations  at  Subsonic,  Transonic,  and 

Supersonic  Speeds."  NACA  Rept.  1307,  1957. 
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(see  fig.  10)  on  the  plane  of  the  horizontal  panels,  an  additional 
upwash  W (yf)  is  obtained  as  follows 

WB(yf'  “ Wo 

yf 

Accordingly,  in  the  application  of  equation  (1-13)  the  value  of  ag 
should  be  taken  as 


as  = 


wi (yf)  + w2 (-yf) 


2V_ 


W 

C 

V 


a 1 Yi 


00  Yc‘ 
s * f 


(1-16) 


The  vertical  fins  produce  a sideforce  which  can  be  calculated  the 
same  way  as  the  normal  force  due  to  the  horizontal  fins.  Let  Vx (-z^) 
be  the  sidewash  velocity  on  the  lower  panel  and  V2(z^)  be  that  on  the 
upper  panel.  The  sidewash  angles  3 for  the  symmetrical  and  unsymmetri- 
cal  parts  of  the  flow  are 


Bs  " 


Vx  (-zf)  + V2(zf) 
2V~ 


0 


u 


V1 (~zf)  - V2 (zf) 
2V 


(1-17) 


The  sideforce  due  to  the  addition  of  the  vertical  panels  is  then 


(C  ) - (C  ) 

. BV  B 


(<) 

i \**/v 


7T 


f (sv  - -a 


/' 


BS(CC<>4  d2f 


(1-18) 


(cCi>.  " 4 


V 


<sv2zf2  - g4)  <sv2  - 2f2) 


2c  2 


(1-19) 


s ^z- 

V f 


The  value  of  f3g  should  be  augmented  for  Beskin  sidewash  as  follows: 


es  - 


V2  (z  ) + Vx  (-zf)  V 


o a* 


2V 


V 


(1-20) 


00  Z* 
S f 


In  determining  the  contributions  of  the  horizontal  and  vertical 
panels  to  the  pitching  moment  My  and  yawing  moment,  Mz  , respectively, 
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we  assume  that  their  respective  loads  are  concentrated  at  the  quarter- 
chord  position  of  the  mean  aerodynamic  chords  of  the  exposed  panels. 

Let  be  the  distance  from  the  store  center  of  moments  to  the  quarter 

chord  of  the  MAC  of  the  horizontal  panels,  and  let  the  corres- 

ponding distance  for  the  vertical  panels  (positive*  for  fins  behind  the 
center  of  moments) . The  contribution  of  the  horizontal  panels  to  the 
pitching-moment  coefficient  is  then 


(C  ) - (C  ) 

. m BH  m B 


M 


q“  SR^R 

S 


(Cz)  - (C.) 

BH  BJ 


~r  (i-2i) 

K 


The  contribution  of  the  vertical  panels  to  the  yawing  moment  coefficient 
is 


M 


<c„>  - (Cn> 

L BV  B- 


^ SR^R 
s 


(1-22) 


1-4.  ROLLING  MOMENT  DUE  TO  EMPENNAGE 

First  consider  the  rolling  moments  due  to  the  horizontal  panels  of 
the  empennage  as  if  no  vertical  panels  existed  . Such  results  will  apply 
to  a store  with  only  two  panels,  the  planar  empennage  case,  because  the 
presence  of  vertical  panels  does  have  the  effect  of  reducing  the  rolling 
moment  due  to  the  horizontal  panels.  This  panel-panel  interference  effect 
is  treated  subsequently  for  the  curciform  case. 

The  rolling  moment  due  to  the  horizontal  panels  is  associated  with 
a of  equation  ( I — 1 ) . The  asymmetric  angle  of  attack  is  equivalent 
to  positive  control  deflection  of  the  right  panel  and  negative  deflection 
of  the  left  panel  to  produce  negative  rolling  moment  in  an  aileron  manner. 
In  order  to  determine  the  rolling  moment,  consider  an  application  of  the 
reverse  flow  theorem,  equation  (1-2) , with 


ad 

ad 


0 on  S_  a = 0 on  S_ 
or  ij 

au  °n  SH  ar  = Yf  °n  s 


P 


d 


P 

r 


Ps 
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We  find  that 


Mx  " 2q«  II  Pdy£dSH  *‘2q«  If  “up=dSH 


(1-23) 


H 


H 


since  the  body  contributes  no  rolling  moment.  Let  (cc  0r  be  the  span 

loading  on  the  body  with  horizontal  panels  in  reverse  flow.  The  local 

angle  of  attack  a = y^  can  be  reproduced  by  a roll  rate  of  p = V 

r r s 

about  the  xf  axis  in  the  clockwise  direction  since  the  local  angle  of 

attack  for  such  a condition  is 


py< 


ar  " 


(1-24) 


Since 


t.e. 


(cci>s  ■ 2 


I 


P5  dx, 


(1-25) 


1 . e . 


we  may  write 


Sh 

I au  <cc£)a 


(1-26) 


The  rolling-moment  coefficient  due  to  the  horizontal  panels  is  thus 


M 


(C.)  = 


1 H ^ SR^R  SR^ 


R 


r 

J a» 


(cct)  dyf 


(1-27) 


where  including  the  roll-damping  term  in  equation  (1-1) 


au  = 


wi (yf ) - w2(-yf)  pyf 


2V„ 


V 


(1-28) 


From  equation  (3-37)  of  reference  16, 
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. . /,  2 -i  2ash  \ / a2 

cc.  =1  + — cos  ) I + — 

1 » \ T S * * a*/  V £ y£ 


Uh2  - yf2) (sh2yf2  - a4) 


s 2V  2 
n yf 


' (y*  ‘ £) 


cosh 


-1 


(y  2 + a2)  (s  2 - a2) 


(y  2 - a2)  (s  2 + a2 )_ 


(1-29) 


In  order  to  scale  (C J down  from  slender-body  theory  to  linear  theory, 

**  H 

we  apply  the  factor  ^ of  equation  (1-12)  to  obtain 


1 Ada  'h 


(C.)  = - “ 

* H (S,  - a)2/ 


R a 


f h 

J au(cci)5  dyf 


(1-30) 


This  formula  will  determine  the  rolling  moment  due  to  the  horizontal 
panels  neglecting  any  interference  effects  due  to  the  vertical  panels, 
the  planar  empennage  case. 

The  interference  effects  due  to  the  vertical  panels  on  the  rolling 
moment  due  to  the  horizontal  panels  can  be  shown  by  the  following 
sketch . 


I 


Let  the  horizontal  panels  have  an  au  distribution  for  which  positive 
pressure  occurs  on  the  upper  surface  of  the  right  panel  and  suction  pres- 
sure on  its  lower  surface.  The  resulting  panel  force  produces  positive 
rolling  moment.  On  the  left  panel  the  signs  are  all  opposite,  but 
positive  rolling  moment  is  generated.  Consider  now  the  right  hand  side 
of  the  upper  panel  which  is  undeflected.  It  will  generally  be  subject 
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to  a positive  pressure  field  associated  with  the  flow  above  the  right 
panel  with  a negative  pressure  field  on  its  left  surface.  The  effect  is 
a negative  rolling  moment  which  opposes  that  of  the  horizontal  panels. 

The  rolling  moment  due  to  the  horizontal  panels  has  thus  been  reduced  by 
the  presence  of  the  undeflected  upper  panels  by  virtue  of  reversed  roll 
developed  by  them  through  interference.  It  may  also  be  argued  that  the 
presence  of  the  vertical  panel  presents  a dam  between  the  first  and  second 
quadrants  which  inhibits  communication  between  these  positive  and  negative 
pressure  regions,  and  that  this  pressure  difference  should  therefore  be 
increased  by  the  presence  of  the  vertical  panels.  While  some  increase  in 
the  rolling  moment  due  to  the  horizontal  panels  may  be  ascribed  to  this 
factor,  it  is  more  than  offset  by  the  reverse  rolling  moment  of  the 
vertical  panels  on  the  basis  of  slender-body  theory. 

In  order  to  determine  the  rolling  moment  of  a cruciform  configu- 
ration with  distribution  au(yf)  on  the  horizontal  panels  and  0u(zf) 
on  the  vertical  panels,  we  write  the  reverse  flow  theorem  for  the  vertical 
and  horizontal  panels  with  the  body  aligned  in  the  free  stream  direction. 


ff  “dPrdSH  + ff  BdPrdSv  * / f “rPddSH  + ff  BrpddSv 


H 


"V 


H 


V 


(1-31) 


Choose  the  direct  and  reverse  flow  conditions  as  follows: 


a,  = a a = yr 
d u r yf 


0 , = 0 0 **  z - 

d u r f 


pd - pd 


pr  = p6 


ff  %P»dSH  + ff  BuP„dSV  " ff  *t PddSH  + / f 


H 


V 


H 


zfPddSv  (1-32) 


V 


In  the  integrals  of  equations  (1-31)  and  (1-32)  the  integrals  have  the 
same  signs  because  an  upper  vertical  panel  set  to  produce  a positive  0^ 
would  have  a positive  au  if  rolled  90  degrees  clockwise  without  a 
change  of  setting.  For  the  reverse  flow  configuration  with  ar  = y^, 
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the  horizontal  surfaces  generate  negative  rolling  moment  (counterclockwise 
looking  forward)  and  the  same  is  true  for  the  vertical  surface  with 
= z^.  Accordingly 


M = 

x 


2q« 


S 


s 


yfpddsH  + 2^ 

s 


t . e. 


a P dx  ^ 
u 0 f 


//  2fPddSV 

sv 


I 


3 dx - dz  _ 
u 6 f f 


(1-33) 

At  this  point,  consider  the  vertical  and  horizontal  panels  to  be  of 
identical  planform  so  that 


M 


(a. 


+ e ) 

u 


t . e. 


P 

6 


dxf  dyf 


(1-34) 


In  the  above  expression  3U  is  evaluated  with  z^ 
span  loading  for  the  reverse  flow  case  above  be 


x'  6 


2 


l.e. 


P0  dXf 


equal  to  y^.  Let  the 


(1-35) 


so  that  the  rolling-moment  coefficient  due  to  both  horizontal  and  vertical 
panels  is 


M 


ici> 


x 


HV 


q«>_SR/R 


Vr 


S 

/ 


(a 


u 


6u> 


(CC£)6 


dY< 


The  quantity  au  is  given  by  equation  (1-28)  and  8U  is  obtained 
equation  (1-17)  by  including  the  roll-damping  term 


0 


u 


V-2f>  - 


V 


2 ( zf  ) 


2V 


pzf 


(1-36) 


from 


(1-37) 
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Equation  (1-36)  is  completely  analogous  to  equation  (1-27)  for  the 

rolling  moment  due  to  the  panels  of  a planar  configuration.  In  this  case 

(ccj  is  the  span  loading  on  the  panels  of  a cruciform  configuration 
6 

rolling  in  a clockwise  direction  facing  forward  at  a rate  equal  to 
p radians  per  second.  The  slender-body  theory  yields2 

, v 8R2  T ^ , -i  /sin  2 o\  ^ ± *_-i  /tan  20\  1 

(cc^e  = ~ [cos  29  tanh  ViiTT7;  ■ cos  tanh  (tiTT^J 

+ ^(kj)  sin  40  - 2kx  cos  Al  K(kx)  Z^^k^j  (1-38) 

where 


R = 


1 

2 


cos 


sh5  (a4  + yf4> 
yf 2 (a4  + sh4) 


(1-39) 


(1-40) 


cos  2y 


2R2 


(1-41) 


kx  = sin  2y 


(1-42) 


K(kx) 


dz 

V 1-kj2  sin  2 if 


complete  elliptic  integral 
of  first  kind 
with  modulus  kx 


(1-43) 


Ax 


(1-44) 


Z(A1,k1)  = Jacobi  zeta  function 


(1-45) 


?Adams,  G.  J.  and  Dugan,  D.  W.:  Theoretical  Damping  in  Roll  and  Rolling 

Moment  Due  to  Differential  Wing  Incidence  for  Slender  Cruciform  wings 
and  Wing-Body  Combinations.:  NACA  Rept.  1088,  1952,  Eq.  (18). 


128 


K(k1 ) ZfA^k^  = K(kx)  ElA^ki)  - E (w/2,k1 ) F(A1,kl) 


(1-46) 


F ( A i , k ^ ) 


r\ 

/ 


dz 


V 1-kj2  sin2  z 


incomplete  elliptic  integral 
of  the  first  kind 


E(ir/2,kl) 


E (Ai , ki ) * 


If/ 2 


2 sin2  z dz  = 


complete  elliptic  integral 
of  second  kind 


(1-48) 


- f 


2 sin2  z dz  = incomplete  elliptic  integral 

of  second  kind 


(1-49) 


In  equation  (1-37)  as  0 — < ► y,  that  is,  — ► a,  the  first  of  the  two 

principal  terms  becomes  indeterminate.  It  is  possible  to  show  that 


Lim^  [cos  20  tanh"3  - cos  2y  tanh"1 

— ► - cos  2y  log  (cos  2y)  (1-50) 

Thus,  the  span  loading  at  the  juncture  between  the  panel  and  body, 
yf  = a,  is 

(cc.)  = — — f-  cos  2y  log  (cos  2y)  +0.5  K(k1)  sin  4o|  (1-51) 

/ 6 7T  L J 


As  equation  (1-36)  stands  it  yields  the  rolling  moment  coefficient 
on  the  basis  of  slender-body  theory.  To  obtain  the  linear  theory  esti- 
mate, we  must  multiply  the  right-hand  side  by  ^ from  equation  (1-12) 
to  obtain 


(C 


HV 


i Vd% 

T (s-a)  2f 


/ 


(aU  + 0U) (CCi)6  dyi 


(1-52) 
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While  the  fraction  of  C , associated  with  a or  3 can  be  considered 

JC  vi 

due  to  asymmetrical  flow  condition  on  the  panels  in  question,  the  effect 
of  any  one  set  of  panels  is  manifest  by  pressures  on  both  sets  of 
panels . 
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APPENDIX  II 


EQUATIONS  OF  MOTION  OF  A RIGID 
BODY  WITH  MASS  AND  INERTIA  ASYMMETRIES 

II-l . INTRODUCTION 

The  purpose  of  this  appendix  is  to  present  the  derivation  of  the 
equations  of  motion  of  a rigid  body  with  mass  and  inertia  asymmetries. 

By  mass  asymmetry  we  mean  that  the  center  of  mass  is  not  at  the  origin 
of  the  coordinate  system  fixed  in  the  body  and  used  in  the  equations  of 
motion.  Generally  the  external  store  will  have  axes  of  geometrical 
symmetry  which  will  be  used  in  determining  the  forces  and  moments  acting 
on  it,  but  the  center  of  mass  in  certain  applications  will  not  lie  on 
these  axes  of  symmetry.  By  inertia  asymmetry,  we  mean  that  the  principal 
axes  of  inertia  do  not  coincide  with  the  geometric  axes  of  symmetry  of 
the  store.  For  these  geometric  axes,  the  tensor  of  inertia  will  thus 
possess  off-diagonal  terms,  the  so-called  products  of  inertia. 

I I -2.  AXIS  SYSTEMS  AND  EULERIAN  ANGLES 

In  the  derivation  of  the  equations  of  motion  consider  first  an 
inertial  set  of  axes,  (£,  tj , £)  , fixed  in  the  aircraft  which  is  in  recti- 
linear flight  at  uniform  velocity.  Then  a system  of  axes  (x,  y,  z) 
fixed  in  the  store  is  defined.  The  x axis  is  forward  along  the  store 
longitudinal  axis,  the  y axis  is  positive  to  the  right  looking  forward, 
and  the  z axis  is  vertical  downward  so  that  a right-handed  system 
exists.  The  tj,  £ system  is  an  airplane  fixed  system  with  £ mea- 

sured forward  along  the  fuselage  longitudinal  axis,  tj  laterally  to  the 
right  and,  £ vertically  downward.  The  origin  is  at  the  fuselage  nose. 

In  order  to  determine  the  orientation  of  the  store  with  respect 
to  the  £,  n,  £ axes,  a system  of  angles  will  be  introduced.  The 
system  of  angles  to  be  used  is  that  shown  in  figure  9.  The  origins  of 
the  inertial  and  store  axis  systems  are  taken  at  the  same  point  to  show 
the  angular  orientation.  The  system  consists  of  three  rotations  in  the 
yaw,  pitch,  roll  sequence.  This  order  is  important.  The  three  angles 
will  be  designated 
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yaw 


y 

9 pitch 
* roll 

II-3 . DIFFERENTIAL  EQUATIONS  FOR  THE  ORIENTATION  ANGLES 

In  specifying  the  position  and  orientation  of  the  store  relative 
to  the  aircraft,  we  will  determine  the  coordinates  £ , rj , and  £ of 
the  origin  of  the  x,  y,  z coordinate  system  and  the  angular  quantities 
Y,  0,  and  ct>.  The  equations  of  motion  will,  however,  be  written  in  terms 
of  the  following  quantities  and  their  time  derivatives  £ , tj,  C,  p,  q,  r 
Here  p,  q,  and  r are  the  components  of  the  store  angular  velocity  about 
the  x,  y,  and  z axes,  respectively,  A set  of  differential  equations 
is  required  for  Y,  0,  and  0 in  terms  of  p,  q,  and  r which  are  to 
be  integrated  along  with  the  equations  of  motion  to  keep  track  of  the 
angular  orientation  of  the  store.  These  equations  are  from  reference  15 


¥ - (q  sin  <t>  + r cos  $)/cos  0 

0 = q cos  # - r sin  $ (II-l) 

<t>  = p + q sin  <t>  tan  0 + r cos  $ tan  0 


With  initial  values  of  the  orientation  parameters  0,  and  $>,  the 

foregoing  equation  can  be  integrated  to  obtain  the  variations  of  these 

quantities  with  time. 

From  the  values  of  Y,  0 , and  <&  it  is  possible  to  determine  the 
orientation  of  the  x,  y,  z axes  with  respect  to  the  tj,  ( in  the 

usual  direction  cosine  way.  That  is 
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where  (ref.  15) 


cos  9 cos  Y 


sin  $ sin  9 cos  ¥ 
-cos  <t>  sin  ¥ 


cos  <J>  sin  9 cos  Y 
+ sin  $ sin  ¥ 


cos  9 sin  ¥ sin  $ sin  9 sin  Y 

+ cos  4>  cos  Y 


cos  $ sin  9 sin  ¥ 
-sin  $ cos  Y 


-sin  9 


sin  <J>  cos  9 


cos  $ cos  9 


These  data  are  also  necessary  in  converting  forces  calculated  in 
the  x,  y,  z system  to  the  t],  £ system. 


I 1-4.  TRANSLATIONAL  EQUATIONS  OF  MOTION 

We  will  now  write  the  equations  of  motion  for  translation  of  the 

origin  of  the  x,  y,  z body  axes.  The  origin  is  not  assumed  to  be  at 

\ 

the  body  center  of  mass.  Consider  the  following  sketch. 
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Let  a be  the  vector  from  the  origin  of  the  inertial  system,  point  P, 
to  the  origin  of  the  body  system,  point  0.  Let  r*  be  the  vector  from 
the  origin  of  the  body  axes  to  the  body  center  of  mass. 

*r‘  = "e x x 4-  ’e y y + "e z z (II-3) 

Let  ~r  be  the  vector  from  the  origin  of  the  inertial  system  to  the 
body  center  of  mass.  Then 

T = "a  + "r ' (II-4) 

If  F is  the  vector  force  acting  on  the  body  the  equation  of  trans- 
lational motion  is 


F = m 


d2r 

dt2 


(II-5) 


This  equation  will  now  be  expressed  in  terms  of  the  force  components 
(F^,  F , Fz)  along  x,  y,  and  z;  the  position  coordinates  (x,  y,  z)  of 
the  body  center  of  mass;  and  the  components  of  the  body  angular  velocity 
a)  given  by 


(D  = e p+e  q+e  r 
x ^ y z 


(H-6) 


To  establish  the  acceleration  vector  equation  (II-4)  is  differentiated 
with  respect  to  time 


dr  _ da  dr  * 

dt  dt  dt 


da 

dt 


+ o>  X r1 


(II-7) 

(11*8) 


Differentiating  again  yields 


d£r 

dt2 


da 

dt2 


da)  _ — • ►.  . — ► dr  * 

dt  X r + “ x St 


(II-9) 
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Since 


do) 

dt 


ex  P + 


e q + 

y 


(11-10) 


we  find 


or 


d2r  = dfa 
dt2  dt2 


— - + ( "e  p + "e  q + e r)  X r 1 + (I  X (la  X 

r\+-2  \ x y z ' V / 


d2a 

dt2 


(** 


p+eyq+ezrj  x r ' + cu 


) x 7‘  + 7(7-7)  - 7'  (7*7) 


d2r  = dfa. 
dt2  dt2 


+ e 


x 


(pq  - r)y  4-  (pr  4-  q)  z - x(q2  4-  r2) 


4*  e 


4-  e 


(pq  4-  r)x  4-  (qr  - p)  z - y(p2  4-  r2) 


(pr  - q)x  + (qr  4-  p)  y - z (p2  + q2) 


(II-ll) 


which  from  equation  (II-5)  is  equal  to 


d2r  = 1 
dt2  m 


Fe  4-  F e 4-Fe 
xx  y y z z_ 


(11-12) 


The  acceleration  of  0 relative  to  P can  be  expressed  in  terms 
of  the  moving  body  coordinate  system  as  follows: 


d2"a 

dt2 


x e 
o x 


ye 

2 o y 


z e 
o z 


+ a)  x 


da 

dt 


(H-13) 


where  — is  the  velocity  of  0 relative  to  P and  may  be  expressed 
in  terms  of  components  along  the  moving  body  coordinate  system 
xQ,  YQ>  and  ZQ*  The  body  angular  velocity  <0  is  given  by  equation 
(II-6) . Therefore 
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d^a  ••  — ► ..  -►  • * x-* 

= xe  + ye  + ze  + (qz  - y r)  e 
(jt2  ox  Joy  oz  ^ o Jfr>  ' 


o x 


+ (x  r - pz  )e  + (py  -xq)e 
o ^ o y o z 


(11-14) 


Combining  equations  (11-11) , (11-12) , and  (11-14)  gives 


F 

X 

m 

£— (q2  + r2)x  + (pq 

- r)  y + (pr  + q)  z 

’]  - (<jzo  - y0r)> 

F 

r . 

• *1 

_ 

m 

[(pq  + r)x  - (p2  + 

r2)  y + (qr  - p)  z J 

\ - (xor  - pzo)  ) 

F 

P • 

z 

m 

[(pr  - q)x  + (qr  + 

p)y  - (p2  + q2)z  J 

1 - (pyo  - xoq)  J 

(H-15) 


The  accelerations  along  the  inertial  axes  are  then 


— ••  —i 

X 

o 

n 

n 

ry. 

• • 

z 

- O J 

where  the  transformation  matrix  is  given  by  equation  (II-2) . 


The  final  equations  of  translational  motion  in  the  inertial  coor- 
dinate system  are  now  written  using  the  following  direction  cosine 
notation  for  the  A matrix  of  equation  II-2 


|"cos 

(£,x) 

COS 

cos 

(n,x) 

cos 

cos 

(C,x) 

cos 

(£,y)  cos  (£,z) 
( T] , y)  cos  (t),z) 

(C,y)  cos  U>z) 
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£ + [y  cos  (£,z)  - z cos  (f,,  y)]p  + [z  cos  (£,x)  - x cos  (£,z)]q 
+ ^ x cos  (£,y)  - y cos  (£,x)]r 

= [-j^-  + x (q2  + r2)  - ypq  - zpr  - qzQ  + ryQ]  cos  (£,x) 


P F . . -j 

+ - xpq  + y (p2  + r2)  - zqr  + pzQ  - rxQJ  cos  (£,y) 


+ - x, 


l m 


xpr  - yqr  + z(p2  + q2)  + qx 


xo  - pyG] 


cos  (£,z) 


(11-16) 


t|  + £y  cos  (tj,z)  - z cos  (r),y)]p  + [z  cos  (q,x)  - x cos  (t),z)J  q 

+ [x  cos  (T),y)  - y cos  (tj,x)]  r 

rF  _ _ . . I 

= + x ( q2  + r2)  - ypq  - zpr  - qzQ  + ryQJ  cos  (q,x) 


+ ["m1  ~ xpq  + y(p2  + r2)  - zqr  + Pz0  " rxo]  cos  ^.y) 


+ r!2  _ 


L m 


xpr  - yqr  + z (p2  + q2)  + qxQ  - pyQ]  cos  (r|,z) 


(11-17) 

C + [y  cos  (C,z)  - z cos  ( C , y ) J P + z cos  (C,x)  - x cos  (C,z)J  q 
+ ^x  cos  (C,y)  - y cos  (£  ,x)]  r 

r F • • , 

= (."rtf  + x(q2  + r2)  “ ypq  “ zpr  “ qzo  + ryoJ  cos  (^’x) 

r f _ . I 

+ - xpq  + y(p2  + r2)  - zqr  + pzQ  - rxQJ  cos  (£, y) 

rFz  _ _ _ • i 

+ - xpr  + yqr  + z (p2  + q2)  + qxo  - pyQJcos  (;,z) 

(11-18) 
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II— 5 . EQUATIONS  OF  ROTATIONAL  MOTION 

The  equations  of  rotational  motion  will  be  written  in  terms  of 
body  axes  because  the  moments  and  products  of  inertia  are  not  functions 
of  time  in  these  axes.  Also,  the  moments  due  to  aerodynamic  forces  are 
calculated  in  this  axis  system  which  is  an  additional  convenience. 
Consider  a general  body  with  axes  x,  y,  z having  its  origin  at  point 
0 and  an  inertial  system  £,  rj,  £ having  its  origin  at  point  P.  Let 
p be  the  position  vector  of  an  element  of  mass  in  the  inertial  system 
and  ~p'  be  the  position  vector  in  the  body  system. 


Let  L be  the  torque  of  the  forces  acting  on  the  body  with  respect 
P 

to  point  P.  Then  the  law  of  the  conservation  of  angular  momentum  in 
inertial  coordinates  is 


(11-19) 


where  the  angular  momentum  vector  is  given  by 


(11-20) 


the  integral  being  taken  over  the  body  mass. 

For  the  store  the  origin  0,  of  the  body  system,  will  generally  be 
offset  from  the  body  center  of  gravity  by  position  vector  r',  and  the 
tensor  of  inertia  for  the  x,  y,  z axes  will  generally  contain  products 
of  inertia.  The  torque  will  be  known  with  respect  to  the  origin  0, 
and  the  angular  velocity  is  given  by  equation  (II-6)  in  terms  of 
components  p,  q,  and  r in  body  axes.  We  will  now  transform  equation 
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(11-19)  to  body  axes  in  terms  of  the  foregoing  quantities.  The  angular 
momentum  becomes 


"h  — J (a  + "p * ) X ("a  + “p 1 ) dm 


Let 


H 


-/■ 


p'x^" dm 


so  that 


H 


"• * / n d™ + 


/?■  - If 


dm  + H 


and 


dH 

E 

dt 


= a x 


f0  dm  + dt  X /dt  dm  + ~dt 


+ fite  dm  x If + /"p’dm  x 


Certain  simplifications  arise  in  the  foregoing  expression  from 
property  of  the  center  of  mass 


so  that 


/ 


(p ' - r ' ) dm  = 0 


/ 


p 1 dm  = r ' m 


and 


dm  = f~?' 1 dm  = m 


Also 


(11-21) 

(11-22) 


(11-23) 

the 

(11-24) 

(H-25) 

(11-26) 


139 


dm  = at?  /'p'  ^ + d±*  /'a  dm 


dar 1 , d2; 

m dta  + m dt2 


m (Sx?-)  + m 


f-*  „ dr’  den  _ d2a 

= mhXlF+dtXr  +dF 


m - T? 

m dt2  “ P 


(11-27) 


Accordingly,  we  find  from  equation  (11-19)  that 


L = a X F + m 
P 


da  (da  d?,N\  dr'  da 

dt  X \dt  + dty  m dt  X dt 


d2!T  dHo 
+ m r X dt^+  “dt" 


which  becomes 


_ d2-‘-  dH 

tp  = ‘a  X F + m r-1  X^Ef+-^- 


(11-28) 


Now  define 


dar 


ro  5 m x 0 + “dF 


(11-29) 


where  Lq  is  the  torque  with  respect  to  the  0 origin  so  that 


L = L - a X F 
o p 


(11-30) 


Equation  (11-29)  for  the  conservation  of  angular  momentum  with  respect 
to  origin  0 has  the  same  form  as  that  with  respect  to  origin  P plus 
an  additional  term. 

With  regard  to  the  three  terms  of  equation  (11-29) , the  torque  Lq 
will  have  three  components  along  the  x,  y,  z axes 


L=eM+eM+eM 
o xx  y y z z 


(II-31) 
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The  acceleration  vector  d2a/dt2  has  components  r),  and  £ which  can 

be  expressed  in  terms  of  components  x , y , and  zq  along  the  x,  y,  z 
axes  as  follows 

••  ..  ..  *• 

= £ cos  (x,£)  + J]  cos  (x , Tj)  + C cos  (x,C) 

yo  * i COS  (y,£)  + T]  cos  (y,T))  + C cos  (y,0  ) (11-32) 

zq  = £ cos  (z,£)  + T|  cos  ( z , T|)  + £ cos  (Z>C)  .> 


The  direction  cosines  in  the  above  equation  are  obtained  from  the  trans- 
pose of  [a] , equation  (II-2) . 

With  the  following  notation 


r ' = 


ex  x + 


e y 

y 1 


+ e 


(11-33) 


we  have 


X 0 = ni?x  [ y(zQ  + pyo  -qxQ)  - z (yQ  + rxQ  - pzQ)  ] 

[•»  • • _ »•  • • 
z(xo  + qzQ  - ryo)  - x(zq  + pyo  - qxo)J 

+ mez  [x(yQ  + rxQ  - pzQ)  - y(xQ  + qzQ  - ryQ)]  (11-34) 


With  regard  to  the  remaining  term  in  equation  (11-29) , we  have 
from  equation  (11-22) 


dt 


dm 


We  first  evaluate  H and  then  differentiate 

o 
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»o  - h- 


x -gt  dm 


= J~~p'  X (m  x p*' ) dm  =J'f'mCp'  • "p ' ) - "p"'  ("p ' • u>)J 


dm 


= I 7o(x2  + y2  + z2)dm  - / (e  x + "e  y + "e  z)  (px  + qy  + rz)dm 

J J x y z 


= e 


£p  (y2  + z2)dm  - q J~  xy  dm  - r J xz  dmj 
2^  £-p  J"  x y dm  + q ^ (x2  + z2)  dm  - r J yz  dmj 

£-p  xz  dm  - q J yz  dm  + r ^ (x2  + y2)  dmj 


+ e2 

The  moments  and  products  of  inertia  are  defined  as  follows: 

I 


xx 


= J*  (y2  + z2)  dm  I'ky  s J xy  dm  ^ 


I = / (x2  + z2)  dm  I = xz  dm 

yy  J xz  J 


[2Z  = / {X2  + y2)  dm  1yz  = / VZ  dm 


(11-35) 


(11-36) 


so  that  equation  (11-35)  becomes 

*o  ex  [^xx  ^*xy  rIxz] 


+ e 


f - pi  + ql  - rl  1 
y L xy  yy  yzJ 


+ e [-pi  - ql  + rl  1 
z L xz  ^ yz  zzj 


(11-37) 


In  body  coordinates  the  moments  and  products  of  inertia  are  not  functions 
of  time  for  a rigid  body.  With  the  help  of  the  relationships 
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de 


df  = ®y  r “ * = “ x ^ 


de 

dt 

dt 


ez  P " ex  r " “ X ey 


e q - e p = a)  x e 
x ^ y ^ z 


(11-38) 


we  find 


dt  = ex  [^>Ixx  " <^Ixy  ~ rIxz  4 r<^^zz  ~ *yy^ 
- (q*  - r»)  I - P(qIJ(2  - rlxy)] 


+ Is  ("-  pi  + ql  - rl  + rp(I  - I ) 
y L xy  ^ yy  yz  ^ xx  zz 


- (r2  - p2)  I - q(rl  - pi  •)] 

* xz  ^ xy  r yz  J 

+ ez  ^*xz  * ^yz  + rIzz  + pq^yy  *xx^ 

- (P2  - q=)Ixy  - r(plyz  - qlxz>]  (11-39) 

Equating  the  components  of  equation  (11-31)  to  the  sum  of  equations 
(11-34)  and  (11-39)  yields  the  three  desired  equations  of  rotational 
motion 
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M 


= ™[y(zo  + pyQ  - q*0)  - z(yc  + rxo  - Pzo)]  + Piyv  - qi„„  - n. 


xx  ^ xy  xz 


M 


+ rq(Izz  - Iyy)  - (q2  - r2)Iyz  “ p(qIxz  " rIxy> 

= m[7(xo  + qzQ  - ryQ)  - x(zo  + pyQ  - qxQ)]  - plxy  + ql^  - rly2 


+ rp(I  - I ) - (r2  - p2) I _ - q (rl  - pi  ) 
r xx  zz  r xz  ^ xy  r yz 


M 


p ••  • • ^ ••  • • 1 • • 

= m x(y  + rx  - pz  ) - y(x  + qz_  - rx  ) - pi  - ql 

l 1 o o ^ o J o ^ o o J r xz  ^ 


+ rl 
yz  zz 


+ pq(I  - I ) - (p2  - q2)I  - r(pl  - ql  ) 
yy  xx  ^ M ' xy  * yz  M xz' 


(11-40) 


The  rotational  equations  of  motion  in  the  inertial  coordinate  system, 
using  equation  (11-32) , are 

m £ y cos  (z,£)  - z cos  (y,^)]^  + m [y  cos  (z,r|)  “ z cos  (y,T])]r| 

+ m [y  COS  (z,C)  - z cos  (y,C)]c  + IXXP  - Ixyq  - Ixzr 

= M - rq (I  - I ) + (q2  - r2)  I + p(ql  - rl  ) 

x M zz  yy'  ^ ' yz  xz  xy' 

- m [y(pyQ  - qxQ)  - z(rxQ  - pzQ)]  (11-41) 

m £z  cos  (x,£)  - x cos  (z,£)Jc  + m £z  cos  (x,r|)  - x cos  (z,q)Jq 

+ m [z  cos  (x,C)  - x cos  (z.C)jc  - IxyP  + Iyyq  - lyzr 

= M - rp(I  - I ) + (r2  - p2)I  + q(rl  - pi  ) 

y xx  zz'  v * ' xz  xy  * yz' 

- m [ z^(qzo  - ryQ)  - x(pyQ  - qxQ)J  (11-42) 
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m [x  cos  (y,£)  - y cos  (x,£)J  £ + m £x  cos  (y,r|)  - y cos  (x,rj)J  r\ 

+ m [x  cos  (y,0  - y cos  (x,C)]  C - IX2P  - Iyzq  + Izzr 

= M - pq  ( I -I  ) + (p2  - q2)I  + r(pl  - ql  ) 
z yy  xx  ^ ^ xy  ^ yz  ^ xz' 

- m [x(rxQ  - pzo)  - y(qzQ  “ r*yo)]  (11-43) 

I 1-6.  OVERALL  SOLUTION  OF  EQUATIONS  OF  MOTION 

Certain  input  data  and  initial  conditions  are  to  be  specified 
before  any  integration  of  the  equations  of  motion  can  be  performed. 

First  the  input  data  must  include  the  body  mass,  m,  and  the  moments  and 
products  of  inertia.  Secondly  the  offset  of  the  origin  of  the  body 
coordinate  system  from  the  center  of  mass  must  be  specified.  We  also 
need  the  initial  body  orientation  as  specified  by  the  three  initial 
angles.  With  regard  to  initial  conditions,  the  initial  position 
(£,  T| , £)  of  the  origin  of  the  body  coordinate  system  and  its 
velocity  (£,  rj , 0 must  be  specified  together  with  the  initial  angular 
velocity  components  p,  q,  and  r.  We  may  then  integrate  equations 
(11-16)  through  (11-18)  and  equations  (11-41)  through  (11-43)  to  determine 
£,  t)?  P,  q>  and  r as  functions  of  time.  Simultaneously  we  must 
integrate  equation  (II-l)  to  keep  track  of  the  body  orientation.  The 
aerodynamic  force  and  moment  components  are  calculated  step-by-step 
during  the  integration  of  the  trajectory. 
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APPENDIX  III 

A REFINED  THEORY  FOR  THE  CALCULATION  OF  NORMAL - 
AND  SIDE-FORCE  DISTRIBUTIONS  ON  BODIES  OF 
REVOLUTION  IN  NONUNIFORM  FLOW 

III-l . INTRODUCTION 

This  appendix  is  concerned  with  a method  for  determining  the  axial 
distributions  of  normal  and  side-force  on  bodies  of  revolution  immersed 
in  nonuniform,  irrotational , incompressible  flow  fields.  This  theory  does 
not  employ  the  slender  body  assumption  nor  the  linear  form  of  the  flow 
tangency  boundary  condition.  It  is,  therefore,  considered  more  refined 
than  the  body  force  calculation  methods  described  in  section  5.2. 

In  the  following  sections  the  general  flow  tangency  boundary  con- 
dition will  be  formulated  and  specialized  for  a body  of  revolution.  It 
will  be  shown  how  the  problem  can  be  represented  as  a superposition  of 
three  different  distributions  of  discrete  singularities  laid  out  along 
the  body  axis.  A boundary  condition  is  applied  at  the  body  surface  for 
each  distribution.  The  flow  tangency  formulations  allow  for  the  deter- 
mination of  the  unknown  singularity  strengths  which  are  necessary  for  the 
computation  of  the  resultant  velocity  on  the  body  surface.  The  pressure 
coefficient  is  then  calculated  from  the  resultant  velocity  and  finally 
the  body  force  distributions  are  expressed  in  terms  of  integrals  of  the 
pressure  coefficient  around  the  surface  of  the  body  of  revolution. 

Effects  of  buoyancy  are  not  accounted  for  in  this  method  but  may  be 
determined  independently  and  their  effects  added  to  the  results  of  this 
appendix . 

III-2 . FLOW  TANGENCY  BOUNDARY  CONDITION  APPLIED  AT  BODY  SURFACE 

The  coordinate  systems  and  part  of  the  nomenclature  used  in  this 
appendix  are  indicated  in  the  following  sketch.  In  the  cylindrical 
coordinate  system  (x,r,0),  the  nonuniform  free-stream  velocity  vector 
is  given  by 

"^(Xjr,#)  = [v(x)  cos  0 + W(x)  sin  0j"eTr 

+ £w  (x)  cos  G + V (x)  sin  + U(x)"ex  (III-l) 
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where  U(x),  V(x),  and  W(x)  are  the  components  of  the  nonuniform  free- 
stream  velocity  in  the  x,  y,  and  z directions,  respectively. 


In  general,  the  body  surface  is  specified  by 


F (x , r , 0)  = 0,  0^x^^c 


Then  the  unit  normal  vector  to  the  body  surface  is  given  by1 


— ► qrad  F 

n = -p3 ■,  „ , 

grad  F 


where 


grad  l'-ff®r  + 7§fee  + f1exj 


(III-2) 


(III-3) 


The  perturbation  flow  field  induced  by  the  body  is  characterized  by  a 
perturbation  velocity  vector  q(x,r,0)  defined  as 


xHarry  F.  Davis,  Introduction  to  Vector  Analysis,  Allyn  and  Bacon,  Inc., 
Boston,  1962,  p.  151. 


147 


(III-4) 


q (x, r , 0) 


ur(x, r,0)  er  + 


u0(x,r,0)e0  + ux(x,r,0)ex 


The  components  of  the  perturbation  velocity  are  u , u , un  in  the 
x,  r,  and  9 directions,  respectively.  The  flow  tangency  boundary  con- 
dition is  satisfied  at  the  body  surface  and  formulated  as  follows 

[ ’woo(xirie)  + ■q(x,r,0)J  • "n  = 0 (III-5) 

where  the  item  inside  the  brackets  represents  the  resultant  velocity 
vector,  \7r.  After  replacing  the  velocity  and  unit  vectors  by  equations 
(III-l) , (III-3)  and  (III-4) , the  result  is 


cos  0 + w(x) 


sin  0 + 


+ V (x) 


sin 


0 + u0 (x, r , 0) j 


ur(x,r,0)j  ^ »-  [w (x)  cos  0 ^ 

? H + [™>i + ux<>'.r.e»]  ° > 

F(x,  r,0)  =0;  0 £ x £ i 

sj 


(III-6) 


Equation  (III-6)  represents  the  non-linear  flow  tangency  condition  for 
any  body  situated  in  a non-uniform  flow.  In  the  treatment  that  follows 
in  this  appendix,  the  only  constraint  will  be  the  condition  of  the  body 
being  a closed  body  of  revolution  represented  by 


F(x,r,0)  = r - R (x)  = 0 or  r = R(x) 
dF  , dF  dR (x)  dF 

SF  7 L’  Fx £r'  'be  = ° 


(III-7) 


The  flow  tangency  boundary  condition  can  therefore  be  rewritten  for  a 
body  of  revolution  as  follows 


ur  (x,  r , 0)  = £u  (x)  + ux(x,r,0)J  - V (x)  cos  9 - W(x)  sin  9 ^ 


r = R (x)  ; 0 £ x £ ic 


Note  that  perturbation  velocity  u^(x,r ,9)  is  not  neglected  as  is  done 
in  linearized  flow  theory. 

The  boundary  conditions  given  by  equation  (III-8)  will  be  split 
in  three  parts  in  such  a way  that  the  problem  of  a body  of  revolution 
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situated  in  a nonuniform  flow  field  may  be  expressed  as  the  superposition 
of  three  separate  problems.  The  three  problems  are: 


(1)  Axisymmetric  flow  past  a body  of  revolution  with  the  flow 

field  characterized  by  velocity  U(x) ; axisymmetric  problem 
with  boundary  condition 

ur  (x,r)  = [u (x)  + ux  (x, r)J 
s s 

r = R(x) ; 0 £ x £ ig 

(2)  Lateral  flow  past  the  same  body  with  the  flow  field  characterized 

by  velocity  V(x) ; sidewash  problem  with  boundary  condition 


u (x, r , 0)  = u (x,r  ,0)  - V(x)  cos 


dx 

'D,  V D,  V 

r = R(x) ; 0 £ x £ 


(III-10) 


(3)  Lateral  flow  past  the  same  body  with  the  flow  field  characterized 
by  velocity  W(x);  upwash  problem  with  boundary  condition 


u (x,  r,0)  = u (x,r,0)  - W(x)  sin  0 


'D,W 


"D,W  dx 

r = R(x)  ; 0 ^ x ^ lc 


(III-ll) 


The  disturbance  flow  fields  induced  by  the  body  in  problems  (1) , (2) 

and  (3)  can  be  represented  by  suitable  distributions  of  point  sources, 
horizontal  doublets  with  their  axes  opposed  to  sidewash  V(x),  vertical 
doublets  with  their  axes  opposed  to  upwash  W(x).  The  perturbation 
velocity  components  urs,  ur  , urD>w,  and  uXs,  uXd>v,  uXd>w  a«d 
U0D  V*  U0D  w are  ;‘-nduced  by  t^ie  sources  (subscript  s)  , horizontal 
doublets  (subscript  D,V)  and  vertical  doublets  (subscript  D,W) . 
According  to  the  superposition  principle, 

■\ 

</>(x,r,0)  = </>s(x,r)  + 0D^v(x,r,0)  + <f>D  (x,  r,  0) 

u ( x , r , 0)  = u (x, r)  + u (x,r,0)  + u (x,r,0) 
s XD,V  D,W 

(Continued  on  next  page) 
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u (x,r,0)  = u (x,r)  + u (x,r,0)  + u (x,r,0)  ) (III-12) 

r s D,V  rD,W 


uQ(x,r,0)  = un  (x,r,0)  + up  (x,r,0) 

D,V  ^D,W 

The  expressions  for  the  velocity  components  are  given  in  the  next  section 
and  then  substituted  in  the  boundary  condition  (III-9)  through  (III-ll) . 


III-2 . PERTURBATION  VELOCITIES  INDUCED  BY  SOURCES  AND  DOUBLETS 
DISCRETE  REPRESENTATIONS  OF  THE  BOUNDARY  CONDITIONS 

In  this  section,  the  three  problems  mentioned  above  will  be  dealt 
with  separately.  The  boundary  conditions  will  be  expressed  first  in 
integral  form  and  then  in  discrete  form  allowing  for  the  determination 
of  the  singularity  strengths  for  each  of  the  three  problems. 

The  axisymmetric  problem  may  be  represented  by  a suitable  distri- 
bution of  point  sources  along  the  axis  of  the  body  of  revolution  in  the 
interval  0 £ x £ i . Let  q = q(x)  be  the  unknown  source  strength 
per  unit  length  of  a continuous  distribution.  The  perturbation  potential 
and  velocity  components  induced  by  the  source  distribution  at  a field- 
point  (x,r,0)  are  given  by:2 


0s(x,r) 


q(£) 


d£ 

[(x  - £)2  + r sj1/2 


*\ 


. , b<t>s  if 

ur  (x’r)  = "Sr"  = 4tt  / 
s J 


q($) 


[(x  - i)s  + r2] 3/2 


d£  > (HI-13) 


30  * f S . 

X ‘ -S’  ' ^ l q'U  '[T*'  - o.- : *.]'■/. 


According  to  the  boundary  condition  (III-9)  associated  with  the  first 
problem,  the  source  strength  must  satisfy  the  following  integral  equation 


2K.  Karamcheti,  Principles  of  Ideal-Fluid  Aerodynamics.  Wiley  and  Sons, 
Inc. , New  York,  1966,  p.  574 
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q<0 


r 

[(X  - $)2  7 r2]3/2 


"N 


U(x)  + 


i 

-L  f 

4t r J 


q($) 


* - 1 


[(X  - {>■  + r*]3'*  dV  d*  > 


dR(x) 


(III-14) 


r = R (x)  ; 0 £ x £ i. 


y"  q(x)  dx  = 0 since  the  body  is  closed-^ 


Explicit  determination  of  q(x)  from  the  above  equations  is  not  readily 
possible.  Instead,  an  approximate  procedure  is  adopted  based  on  an  axial 
distribution  of  a finite  number  of  discrete  point  sources.  The  no-flow 
boundary  condition  is  applied  at  a finite  set  of  control  points  on  the 
body  surface.  The  strength  of  the  discrete  point  source  is  related  to 
the  continuous  source  strength  distribution  as  follows 


Hs  = q(?)  A£ 

i 


( III-15) 


Then  the  perturbation  velocities  induced  at  a field  point  (x,r,0)  by  a 
finite  set  of  point  sources  are  obtained  by  summing  their  individual 
effects.  For  n discrete  sources  the  results  are 


n 


vx,rl  s [<x . ,,)i * 

1=1  11  J 


n 


UX  (x’r)  “ E 4 


X - X. 


. Si  f (x  - x . ) 2 + r2] 
1=1  L 1 J 


3/2 


(III-16) 


Substituting  these  expressions  into  equation  (III-9)  and  applying  the 

xk>  rk  = R(xk)» 

k = l,2,...n-2]  located  on  the  surface  of  the  body  of  revolution  yields 
the  first  condition  for  the  axisymmetric  problem. 


boundary  condition  at  n - 2 control  points 
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— t 

4tt  Z-. 


tzai  k - v2  * R2<>'k']va 


R(xk)  - (xk  - X.) 


dR(xk) 

dx 


dR(xk) 


' “<V  dx 


(III-17) 


k = 1,2,... n-2 


The  second  condition  is  given  by  the  fact  that  the  body  is  closed. 


n 


1=1 


>V  = 0 

1 


( II 1-18) 


To  allow  for  the  existence  of  a stagnation  point  at  the  nose  of  the  body, 
the  component  of  the  total  velocity  in  the  x-direction  must  equal  zero 
at  x = 0.  Since  that  component  is  given  by  the  sum  of  U(x  = 0)  and 
uXg(x  =0,  r = 0) , equation  (III-16) , the  third  condition  is 


n \i 


^ - u<*  - o) 


( III-19) 


1=1 


The  above  three  equations  constitute  a set  of  n equations  in  n unknown 
source  strengths  from  which  the  |±s^  values  of  the  discrete  point  sources 
can  be  determined.  Body  radius  R and  the  longitudinal  component  U of 
the  nonuniform  free  stream  V must  be  specified  at  stations  x,  along 
the  body  axis,  assuming  there  is  negligible  variation  of  U in  the 
lateral  directions . 

The  second  or  sidewash  problem  is  represented  by  a distribution  of 
horizontal  doublets  along  the  axis  of  the  body  with  the  axis  of  the  doub- 
lets opposing  the  direction  of  the  nonuniform  sidewash  V(x) . Let 
Hy  = ^(x)  be  the  unknown  strength  per  unit  length  of  the  continuous 
doublet  distribution.  The  perturbation  potential  and  velocity  components 
induced  by  the  doublet  distribution  may  be  obtained  from  those  given  in 
reference  cited  immediately  above  through  a simple  coordinate  trans- 
formation on  the  angle  0 . The  results  are 
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r s 

/ II  ( u \ 

r cos  6 

j V^} 

0 

[(x  - i)2  + r2] 

3/2 


as 


~\ 
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Substitution  of  the  above  expressions  for  the  perturbation  velocities  in 
the  boundary  condition  (III-10)  would  result  in  an  integral  equation 
from  which  it  is  not  readily  possible  to  determine  the  unknown,  continuous 
doublet  distribution  ^(x)  . Instead,  an  approximate  scheme  can  be 
devised  based  on  a distribution  of  discrete,  horizontal  doublets  on  the 
body  axis.  The  strength  M-d  • V of  such  a doublet  is  related  to  the 

^ I * 

continuous  doublet  distribution  for  the  sidewash  problem  according  to 

uD.,v  = PvU)  M (HI-21) 

The  perturbation  velocities  induced  at  a field  point  (x,r,0)  by  a finite 
set  of  horizontal  doublets  are  obtained  by  summing  their  individual 
effects.  For  n discrete,  horizontal  doublets  the  results  are 
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Insertion  of  these  expressions  into  equation  (III-10)  and  applying  the 
boundary  condition  at  n points  [x^>  = R(x^)  , k = 1,2, ...nj 

distributed  on  the  surface  of  the  body  of  revolution  yields 
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°lxkli!s;  k = 1,2,  . . .n 

u 

The  above  expression  represents  a set  of  n equations  in  n unknown 
doublet  strengths  from  which  the  Ud^,V  values  of  the  discrete  doublets 
can  be  determined.  Body  radius  R and  the  sidewash  component  V of 
the  nonuniform  free  stream  must  be  specified  at  stations  along 

the  body  axis,  assuming  there  is  negligible  variation  in  the  value  of 
V from  the  body  axis  to  its  surface. 
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The  third  or  upwash  problem  is  represented  by  a distribution  of  verti- 
cal doublets  along  the  body  axis  with  the  doublet  axis  opposing  the 
direction  of  the  nonuniform  upwash,  W(x).  Let  = 

known  strength  per  unit  length  of  the  continuous  doublet  distribution. 

The  perturbation  potential  and  velocity  components  induced  by  this  con- 
tinuous doublet  distribution  may  be  obtained  directly  from  the  reference 
cited  previously. 
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Substitution  of  the  above  expressions  for  the  perturbation  velocities  in 
the  boundary  condition  (III-ll)  would  again  result  in  an  integral  equation 
from  which  the  unknown,  continuous  doublet  distribution  M^(x)  cannot 
be  readily  determined.  As  for  the  case  discussed  above,  this  problem 
can  also  be  solved  by  means  of  an  approximate  method  based  on  a distri- 
bution of  discrete,  vertical  doublets  on  the  body  axis.  The  strength 

UD • W °f  such  a doublet  is  related  to  the  continuous  doublet  distribution 
1 > 

for  the  upwash  problem  as  follows 
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(III-25) 


155 


The  perturbation  velocities  induced  at  a field  point  (x,r,0)  by  a finite 
set  of  vertical  doublets  are  obtained  by  adding  their  individual  effects. 
For  n discrete,  vertical  doublets  the  results  are 
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Inserting  these  expressions  into  equation  (III-ll)  and  applying  the 
boundary  condition  at  n points  x^,  r^  = R(x^)  , k = l,2,...nj  distri- 
buted on  the  surface  of  the  body  of  revolution  yields  the  following  result 
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The  above  expression  represents  a set  of  n equations  in  n unknown 
doublet  strengths  from  which  the  M-Dj^W  values  of  the  discrete,  vertical 
doublets  can  be  calculated.  Body  radius  R and  the  upwash  component 
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W must  be  specified  at  stations  x^  along  the  body  axis,  assuming 
there  is  negligible  variation  in  the  value  of  W from  the  body  axis  to 
its  surface. 

Once  the  strengths  of  the  discrete  sources,  horizontal  and  vertical 
doublets  have  been  determined,  the  resultant  velocity  and  the  pressure 
coefficient  can  be  calculated  on  the  body  surface  by  an  application  of 
the  superposition  principles  stated  in  equation  (III-12)  as  shown  in  the 
next  section. 


III-4.  PRESSURE  COEFFICIENT  ON  BODY  SURFACE  AND  NORMAL-  AND  SIDE- FORCE 
DISTRIBUTIONS  ALONG  BODY  AXIS 

The  pressure  at  any  point  on  the  body  surface  is  related  to  the 
resultant  velocity  through  the  Bernoulli  equation 
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(III-28) 


The  reference  velocity  is  taken  as  the  velocity  the  body  (store)  sees 
at  infinity,  denoted  Vco  . The  nonuniform  free  stream  velocity  "v 
characterizes  the  local  flow  field  in  which  the  body  is  immersed,  and 
is  given  by  equation  (III-l) . The  resultant  velocity  is  obtained 

by  adding  the  components  of  the  free  stream  to  the  appropriate  pertur- 
bation velocities  induced  by  the  sources  and  the  horizontal  and  vertical 
doublets  in  accordance  with  the  superposition  principle  (III-12) . 
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The  magnitude,  VR,  of  the  resultant  velocity  is  consequently  given  by 
the  square  root  of  the  sum  of  the  squares  of  the  coefficients  in  the 
above  equation. 

The  rate  of  change  of  normal  and  side  force  with  x can  be  obtained 
from  reference  1,  equation  (85).  In  conjunction  with  equation  (III-28)  , 
the  results  are 
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After  replacing  the  perturbation  velocities  in  equation  (III-29) 
by  their  discrete  representations  as  given  by  equations  (III-16) , (III-22) 

and  (III-26)  and  substituting  the  result  in  equation  (III-30) , the  surviving 
terms  in  the  integrands  are  those  involving  sin2  6 and  cos2  6 . 

The  integrals  therefore  reduce  to 
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,r,9)  + 2u  (x, r) u (x,r,e) 
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Substitutions  from  equations  (III-16) , (III-22)  and  (III-26)  yield  the 
final  results  for  the  axial  distributions  of  the  normal  and  side-force 
coefficients  acting  on  a body  of  revolution  with  radius  R(x)  and 
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inserted  in  a nonuniform  flow  field  characterized  by  U(x)  , V(x)  and 
W(x)  . 
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In  the  equations  above,  it  has  been  assumed  that  the  number  of  discrete 
sources,  vertical  doublets,  and  horizontal  doublets  are  the  same.  This 
condition  can  be  relaxed  as  long  as  the  number  of  singularities  is 
compatible  with  the  number  of  control  points  in  each  of  the  three  pro- 
blems discussed  in  section  (III-2) . 
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